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CHAPTER I 


INTRODUCTIO N 


1 . 1 General Statement of the Problem 

Turbulent flows containing a dispersed phase, e.g., particles, 
drops, or bubbles, are often encountered during the design of heat 
transfer, mass transfer, and combustion equipment. For example, in 
gas turbine combustors, liquid fuel is generally sprayed into a 
turbulent swirling recirculating flowfield where the liquid breaks up 
into droplets, evaporates, and reacts. Because of the high cost of 
designing combustors, there is substantial interest in developing 
computer models of this very complicated two-phase flow. A better 
understanding of the fundamental processes involved in these flows can 
be used to improve current computer models and subsequently enhance 
the ability to accurately predict flow properties. 

While understanding the entire spray combustion problem is the 
eventual goal, a simpler problem is considered here. The present 
study is limited to the investigation of particle-laden weakly 
swirling free jets. For this case, the only interaction between the 
two phases is the exchange of momentum. The arrangement considered 
consisted of axisymmetric weakly swirling air jets containing 
particles, which are injected downward into stagnant air. The swirl 
number was limited to 0.4 or less so that a recirculation zone is not 
formed and the boundary-layer approximations can be made. Major 
features of the flow that are of interest during the present 
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investigation include: the velocity distribution of both phases, the 

distribution of void fraction, and particle concentrations. The 
continuous phase is turbulent; therefore, local velocity fluctuations 
and Reynolds stresses are important parameters since they control the 
rate of spread of the continuous phase and the entrainment of ambient 
air into the jet. The particles interact with the turbulent 
continuous phase and exhibit random processes as a result; therefore, 
both the mean and fluctuating properties of the particle phase are of 
Interest. Turbulent dispersion of particles, i.e., the motion of 
particles due to their interaction with turbulent eddies, is important 
since it strongly influences the spread of the particles. Finally, 
particle concentrations are valuable since they reflect the mixing 
properties of the flow. 

The present study attempts to resolve some of these features 
emphasizing new measurements in weakly swirling particle-laden flows. 
Predictions are also considered, both to help interpretation of the 
measurements and to initiate evaluation of methods to analyze flow 
properties. The following section will briefly describe previous 
studies of two-phase flows which are related to the present study. 

1 .2 Related Studies 

1.2.1 Single-Phase Swirling Jets 

Single-phase swirling jets have been studied by a large number of 
investigators. Earlier experimental studies of swirling jets have 
been reported in Refs. 1 - 13. Theoretical analysis of some aspects 



of swirling flows is described by Murthy (14)*. The present 
discussion will primarily concentrate on previous studies of 
unconfined swirling jets. 

In the experimental studies reported in Refs. 1 - 12, swirl has 
been generated by various methods. Rose (1) and Pratt and Keffer (2) 
generated a swirling jet by rotating a round tube which was passing a 
fully developed flow. Mathur and Maccalum (3), Kerr and Fraser (4), 
and Sislian and Cusworth (5) used swirlers to generate the tangential 
component of velocity for their swirling flows. Tangential air 
injection into an axial flow has also been frequently used (6 - 11). 
Gouldin, et al . (12) utilized both swirlers and tangential air 
injection. For the earlier experimental studies (3, 4, 6 and 9), 
pitot probes were used to measure velocities; therefore, only mean 
velocities were reported. Hot-wire anemometry has also been used to 
measure both mean and fluctuating velocities (1, 2, 7, 8, and 11). 
Later studies used nonintrusive laser-Doppler anemometry (LDA) to 
measure mean and fluctuating quantities (5 and 12). 

A measure of the amount of swirl in a flow is given by the swirl 
number. The swirl number is the axial flux of angular momentum 
divided by the axial flux of axial momentum multiplied by a 
representative length (15). For a constant density free jet in 
stagnant surroundings, neglecting effects of fluctuating velocities, 
the swirl number can be calculated from the following expression: 

‘Numbers in parenthesis denote references. 
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The swirl number, S, is a conserved property of the flow if the 
ambient conditions are stagnant, i.e.. if the ambient turbulence 


intensity and angular velocity are zero and there is no imposed 
longitudinal pressure gradient. As discussed in Ref. 15, the swirl 
number is an important parameter of a swirling flow. For swirl 
numbers less than about 0.6, the adverse pressure gradient caused by 
the tangential velocity decay is not strong enough to cause a reversal 
of the axial velocity. Such flows are called weakly swirling flows 
and have been analyzed using the boundary-layer approximations (10 and 
11). As the swirl number increases, the rate of growth, the rate of 
entrainment of ambient fluid, and the rate of axial velocity decay all 
increase. At swirl numbers above about 0.6, strong axial and radial 
pressure gradients are present which cause a recirculation zone to be 
formed along the axis. As discussed in Ref. 15, these flows are 
governed by elliptic partial differential equations. 

Lee (16) reports an early theoretical analysis of unconfined, 
asymmetric, turbulent swirling jets, based on similarity of axial and 
swirling velocity profiles, as well as lateral entrainment. The 
comparison between predictions and the measurements of Rose (1), for 
mean values of axial and angular velocity at an axial distance of 
three pipe diameters downstream of the jet exit, were reasonably 
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good. Several computer models of turbulent swirling jets have been 
reported (10, 11, and 17 - 22). A study, reported by Siddhartha (10), 
used a parabolic marching procedure with a turbulence model based on 
the Prandtl mixing-layer hypothesis to predict both free and confined 
weakly swirling jets. Lilley (18) used an identical parabolic code 
with a mixing-length turbulence model that was empirically modified to 
account for anisotropy. Later studies have employed more advanced 
turbulence models. Koosinlin and Lockwood (22) used an algebraic- 
stress turbulence model to predict swirling boundary-layer flows. 

Morse (11) predicted the structure of weakly swirling turbulent free 
jets with a full Reynol ds-stress turbulence model. His predictions, 
using the full Reynol ds-stress model, did not show any better 
agreement with experimental data than much simpler mixing-length 
formulations. This Reynolds-stress closure has recently been 
modified, however, by Gibson and Younis (23) to improve its 
performance in swirling jets. Elliptic type calculations of swirling 
flows, using the two-equation k-e turbulence model, can be found in 
Refs . 1 7 and 19 - 21 . 

Swirling flows are very important in combustion applications. 
Studies of combusting swirling flows can be found in Refs. 13 and 24 
26. 

1.2.2 Particle-Laden Flows 

Previous research on particle-laden flows is quite extensive. A 
comprehensive discussion early work in this area can be found in Soo 


5 



(27). More recent reviews of par t i c I e- 1 aden flows can be found in 
Refs. 28 - 30. Since experimental studies of two-phase flows are very 
numerous, the present discussion will be mainly limited to recent 
studies of particle- or droplet-laden jets. Yuu et al. (31) studied 
particle-laden jets containing fly-ash particles (15 to 20 pm in 
diameter) injected into stagnant air from a nozzle designed to produce 
a uniform outlet velocity. Gas-phase mean velocity and particle 
concentration measurements were reported. McComb and Salih (32 and 
33) measured particle concentrations of 2.3 and 5.7 pm diameter 
particles injected into stagnant air, using 1 aser-Doppl er anemometry 
(LDA) techniques. Popper et al . (34) used LDA techniques to measure 
velocities of oil droplets whose diameters were estimated to be less 
than 50 pm: the jets of oil droplets and air were injected into still 

air from a nozzle designed to produce a uniform velocity at the exit; 
the mass-loading ratio of droplets to air for their flow was limited 
to an extremely low value of 0.001. Levy and Lockwood (35) measured 
mean and fluctuating velocities of both the gas phase and particles, 
using LDA, for a round jet discharging into still air: the injection 

pipe was 40 diameters long in order to obtain approximately fully 
developed turbulent flow at the exit; the particles used were 
relatively large sand particles, in the range of 215 to 1060 pm in 
diameter; while mass-loading ratios were fairly high and ranged from 
1.1 to 3.5. Modarress et al . (36 and 37) used a two-color LDA to 
measure mean and fluctuating velocities of both the particle and gas 



phases: their experimental configuration consisted of a round jet 

discharging into a low velocity co-flowing airstream within a duct, 
with particles consisting of glass beads, 50 and 200 pm in diameter, 
and mass-loading ratios of 0.32 to 0.85. Shuen et al. (38 - 40), 

Shuen (41) and Zhang et al . (42) reported measurements of a 
particle-laden round jet discharging into a stagnant environment using 
three sizes of sand particles (79, 119, and 207 pm in diameter) and 
various loading ratios, while employing LDA to measure the velocities 
of both phases. Solomon et al. (43 and 44) reported measurements in a 
nonevaporating spray of vacuum-pump oil from an air-atomizing 
injector, using a double-flash photographic technique to measure 
droplet velocities and LDA to measure mean and fluctuating velocities 
of the continuous phase. Evaporating sprays were studied by Solomon 
et al . (45), while measurements of combusting monodisperse droplets 
were reported by Shuen et al. (46 and 47): LDA was used to measure 

velocities of both the droplets and the continuous phase in both these 
studies. 

The above experimental studies show that particles or drops exert 
an influence on the continuous phase, decreasing the spreading rate of 
the jet and the centerline velocity decay of the continuous phase. 

The size of the particles, as well as the loading ratio, had an effect 
on the interaction between the turbulent continuous phase and the 
particles . 
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Previous studies of particle-laden swirling flows are relatively 
scarce, e.g., no experimental studies could be found in the 
literature. A few numerical studies have been reported (48 to 51). 
During and Suo (48) obtained solutions for particle trajectories in a 
free-vortex swirling flow. Domingos and Roriz (49) predicted the 
trajectories of evaporating or burning droplets in known gas 
flowfields. Seleznev and Tsvigan (50) performed numerical 
computations for a swirling gas with condensed droplets in an 
expanding channel, treating the gas and droplets as interpenetrating 
media. Finally, Hamed (51) reports particle trajectory calculations 
in a flow field with swirling vanes. Effects of particle impacts 
against the swirler vane surfaces were considered, however, 
interactions between particle trajectories and the gas phase were 
ignored. More details concerning computer models of two-phase flow 
will be discussed in the next section. 

1.2.3 Analysis of Particle-Laden Flo w 

Various computer models have been developed for two-phase flows. 
One-dimensional and lumped-parameter models are discussed in an 
earlier review by Faeth (52) and will not be considered here. 
Subsequent reviews of recent two-phase flow models, applicable to 
sprays, have also been presented by Faeth (53 and 54). 

1 . 2 . 3 . 1 . LHF model s . - There are two major types of two-phase 
flow models. The first involves use of the local ly-homogeneous-f low 
(LHF) approximation where the drops or particles and the gas phase are 
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treated essentially as a single phase. Interphase transport rates are 
assumed infinitely fast under the LHF approximation; therefore, both 
phases have the same temperature and velocity and are also in phase 
equilibrium at every point in the flowfield. The LHF model 
corresponds to a single-phase variable density flow. The LHF 
approximation is strictly accurate only at the limit of infinitely 
small particles. 

LHF models have been extensively used to predict two-phase flows 
because of their relative simplicity and the modest information 
required to specify initial conditions, e.g., detailed specifications 
of initial conditions for the particles are not required. 

LHF analysis of a variety of two-phase flows can be found in 
Refs. 38-47 and 52-56. In general, effects of finite interphase 
transport rates usually resulted in LHF model predictions that 
overestimated the rate of development of the flow. LHF model 
predictions were only satisfactory for flows containing small tracer 
particles, where characteristic response times of the particles were 
small in comparison to all characteristic response times of the 
continuous phase. 

1 .2.3.1 . SF models . - The second type of model considers finite 
transport rates between the phases and is termed a separated-flow 
model. Many separated-flow (SF) models have been proposed. A few are 
described in Refs. 38 - 49 and 51 - 54, and 57 - 73. 
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According to a recent review by Crowe (57), separated-f low models 
of dilute par t i c 1 e- 1 aden gas flows can be divided into two groups, 
two-fluid and lagrangian models. Two-fluid models regard the 
conveying and particulate phases as two interactive fluids similar to 
the two species in a binary mixture. A major disadvantage of this 
approach is that in order to consider more than one particle size, 
each size category must be treated as a separate fluid. This can 
demand a considerable amount of computer storage. An example of the 
two-fluid approach is given by Melville et al. (58). The Lagrangian 
approach is based on the "Parti cle-Source-In-Cel 1 " (PSIC) technique 
reported by Sharma and Crowe (59). The PSIC approach is based on 
treating the particles as sources of mass, momentum, and energy to the 
gaseous phase. Particle trajectories are calculated in conjunction 
with a Eulerian continuous-phase solution. 

An advantage of Lagrangian trajectory calculations is the absence 
of numerical diffusion of the particles. If the Lagrangian 
calculation is carried out using mean values of the continuous phase, 
the particles follow deterministic trajectories similar to the 
behavior of particles in a laminar flow. This type of analysis is 
termed a deterministic separated flow (DSF) model. Examples of DSF 
models are given in Refs. 38 - 47 and 53 - 54. In general, DSF models 
tend to underestimate the rate of dispersion of the particle phase. 

In addition, no information regarding the statistical properties of 
the particle phase is provided by a DSF model. 
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Various methods have been considered in order to treat the 
dispersion of particles by the turbulent fluctuations of the 
continuous phase. One method, first proposed by Jurewicz (60), is to 
use an effective diffusion force that is dependent on the particle 
concentration gradient in the particle motion equation. This requires 
an effective diffusion coefficient for which no reliable information 
is currently available. Another technique used to model the turbulent 
dispersion of particles is based on Monte Carlo methods and is termed 
the stochastic separated flow (SSF) method. In the SSF approach, the 
turbulent gas-phase flow field is modeled as a steady flow with 
superimposed random velocity fluctuations. Particle trajectories are 
then computed by randomly sampling the turbulent properties of the 
gas-phase flow field. This procedure eliminates the need for an 
effective diffusion coefficient, however many particles must be 
considered in order to provide a valid statistical sample. 

A number of methods have been employed to model the mean flow and 
turbulence properties of the continuous phase for SSF models. Some 
applications of the SSF approach used simplified descriptions or 
empirical correlations of turbulence properties (61 - 63). More 
comprehensive treatments have utilized a turbulence model to provide 
the instantaneous properties of the continuous phase. Peskin and Kau 
(64) simulated a particle-laden flow in a rectangular channel using a 
large-eddy simulation (LES) model. The LES turbulence model involves 
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simulation of the three-dimensional, time-dependent, Navier-Stokes 
equations which requires substantial computer time. The well-known 
two-equation k-e turbulence model, developed at Imperial College 
(74), has been used most often (38 - 47, 53, 54, 65 - 67). This is 
due to the lack of versatility of algebraic models, the substantial 
computational requirements of LES, and the relatively early state of 
development of algebraic stress and full Reynolds-stress turbulence 
models. A coherent body of research on two-phase flows is reported in 
Refs. 37 - 47. SSF model predictions based on the k-e turbulence 
model and the turbulent dispersion approach first proposed by Gosman 
and Ioannides (66) are presented for a variety of two-phase flows. 

The present work extends the methods reported in Refs. 37 - 47 to 
weakly swirl ing jets . 

1.3. Specific Problem Statement 

Information on the structure of two-phase, turbulent, weakly 
swirling jets is needed to improve our understanding and advance 
methods of analysis of two-phase flows. 

The previous review of the literature has shown that significant 
progress has recently been made toward developing methods of 
estimating the properties of two-phase flows. However, the rate of 
development of two-phase flow models has outpaced available 
experimental data. In particular, no data are currently available for 
two-phase flows with swirl. Since swirling flows have many important 
practical applications, this is an important area of research. 



The present study seeks to provide new information concerning 
swirling particle-laden flows. Measurements of mean and fluctuating 
properties of both the continuous and dispersed phases were undertaken 
in turbulent, weakly swirling jets. Predictions were also undertaken, 
based on locally homogeneous flow, deterministic separated flow and 
stochastic separated flow models that are typical of their current 
state of development. The predictions were used to help interpret the 
measurements and to initiate evaluation ot methods to analyze swirling 
particle-laden flows. Specific objectives of the study were as 
fol lows : 

1. Measure the structure of particle-laden weakly swirling jets 
injected into stagnant air. Measured properties include: mean and 
fluctuating velocities of both phases, and the distribution of 
particles in the flow field. 

2. Modify existing locally homogeneous flow, deterministic 
separated flow and stochastic separated flow models of dilute 
particle-laden flows so that they can be applied to weakly swirling 
flows . 

3. Compare model predictions with experimental data obtained 
during the present study and use these results to interpret 
measurements and to guide further model development. 

Measurements were limited to relatively monodi sperse , dilute, 
particle-laden free jets. This allowed the use of nonintrusive 
laser-based methods for velocity measurements of both phases, since 
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spurious scattering of the laser beams was modest and particle size 
was known. Boundary conditions for the experiment were well defined. 
Experiments were conducted by injecting an air jet containing solid 
glass spheres downward from a long tube into stagnant air at ambient 
temperature and pressure. Swirl was imparted to the jet by tangential 
air injection and was varied in order to determine its effect. 

Initial conditions for both phases were measured as completely as 
possible, to facilitate calculations of flow properties. 


CHAPTER II 


EXPERIME N TAL ME THODS 


2.1 Introduction 

Measurements were obtained in both single-phase and 
particle-laden jets. The flows examined during the present study have 
properties as follows: 

(1) Three isothermal single-phase jets (with swirl numbers of 0, 
0.19, and 0.33). 

(2) Three particle-laden jets (with a particle SMD of 39 pm, 
loading ratio of 0.2, and swirl numbers of 0, 0.16, and 0.3). 

A variety of measurements were undertaken in order to study the 
flows. Axial and radial profiles of mean and fluctuating velocities 
for both phases and mean particle mass fluxes were measured. Initial 
conditions were measured at x/d = 0.5. Additional measurements 
involved particle size distributions. 

2.2 Experimental Apparatus 

2.2.1 General Arrangement 

A schematic drawing of the experimental arrangement is 
illustrated in figure 2.1. The flow consists of a turbulent, weakly 
swirling particle-laden jet discharging vertically downward into a 
stagnant environment. Swirl numbers of the jets were less than 0.4, 
calculated from equation (1.1), in order to avoid the appearance of a 
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Figure 2.1. - Experimental arrangement. 
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recirculation zone (15). As a result, the flows could be treated 
using the boundary-layer approximations. The boundary-layer 
approximations have previously been applied successfully to these 
flows (10, 11, 15, 16, 18 and 22), even though initial rates of radial 
spread are relatively high. This arrangement provides an experiment 
with a relatively simple geometry and well-defined boundary 
conditions, which is also computationally tractable. 

The single-phase and particle-laden jets discharged from a round 
tube (19 mm inside diameter) that was 1900 mm long, in order to 
provide a fully developed turbulent flow at the discharge plane. The 
injection tube was mounted inside a cage that was 1.8 m square at the 
base and 2.4 m high. The cage was completely enclosed with one layer 
of 16 mesh screen (wire diameter 0.41 mm, open area 55.4 percent) to 
protect the jet from room disturbances. The jet was exhausted to the 
outside using a roof-mounted, variable-speed blower whose inlet was 
roughly 1 m below the measuring plane. The exhaust line (250 mm in 
diameter) began 350 mm below a 60 mesh screen (1.2 x 1.2 m screen with 
wire diameter of 0.28 mm and an open area of 11.7 percent) which acted 
as the inlet to the exhaust system. Velocity measurements at the 
measurement plane, with the exhaust system operating and no jet flow, 
indicated an induced velocity of less than 0.2 m/s from the exhaust 
blower; therefore, it was concluded that the exhaust system had a 
negligible effect on measurements within the jet flows. 
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Most measurements involved optical diagnostics which were rigidly 
mounted; therefore, the entire cage and injector were traversed to 
measure properties at various points in the flow. The cage was 
mounted on two sets of linear bearings in order to provide movement 
along two axes in the horizontal plane. Positioning was accomplished 
using two, 1220 mm long. Unislide assemblies (B6048P20J), driven by 
stepper motors (Slo-Syn M093). Positioning accuracy was estimated to 
be ±0 . 1 mm. To measure properties at various axial locations, the 
tube from which the jet discharged was traversed vertically using a 
1.3 m long Unislide (B6051P20J) driven by a stepper motor (Slo-Syn 
M093). Positioning accuracy in the vertical direction was also 
estimated to be ±0.1 mm. All three stepper motors were controlled by 
a microcomputer driver (Velmex 8300). 

The particles used in the experiments were solid glass spheres 
manufactured by N.T. Ruddock Company. They had a density of 
2500 kg/m 3 and a refractive index of 1.51. The particles were 
manufactured for use in this study and were used as received. The 
particle size distribution for a sample of 1567 particles is 
illustrated in figure 2.2. Particle size was measured using a 
Cambridge Instruments Quantimet 900 image analyzer. The calculated 
SMO for this distribution was 39 pm with a standard deviation of 
15 pm. The average mean diameter of this distribution was calculated 
to be 30 pm. 
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2.2. 2 Flow system 

A sketch of the flow system is illustrated in figure 2.3. 
Compressed air, supplied at approximately 860 kPa, was passed through 
two coalescing filters (Balston 62A-1), mounted in series, to remove 
contaminants. The flow was then separated into three streams that 
were used for the main airflow, the swirl airflow, and a flow which 
was seeded for laser velocimetry (LV) measurements, the flow rates of 
the main and swirl air streams were measured using calibrated, round, 
critical-flow orifices. Various combinations of four orifices (2.286, 
1.27, 1.092 and 0.762 mm in diameter) were used to provide the 
required flowrates. The air flow rate was adjusted by varying the 
pressure on the upstream side of each orifice, using pressure 
regulators (Conoflow Model H40-Xt-HXA). Upstream pressures were 
measured using calibrated strain gage pressure transducers (MBIS 
C-64952-D, 0-689.5 kPa). The static pressure drop across each orifice 
was measured using calibrated strain gage differential pressure 
transducers (Bell and Howell, 4-351-0210, ±689.5 kPa) and checked to 
insure that each orifice was choked at all operating conditions. The 
air temperature upstream of each orifice was measured using 3.2 mm 
diameter closed-end Chroma 1 -A lumel thermocouples, inserted into the 
center of the tube and connected to a Doric Scientific Series 400 
readout . 

The main air flow was seeded with aluminum oxide particles, 
nominally 1 pm in diameter, (Baikalox 1.0 CR) to allow LV measurements 
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Figure 2.3. - Sketch of the flow system. 
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of continuous phase velocities. For the single-phase jets, the 
seeding particles were mixed into the airflow using a f I u i d i zed- bed 
particle generator (TSI Model 3400). For the particle-laden flows, a 
reverse cyclone seeder, described by Glass and Kennedy (75), was used 
in order to provide increased seeding levels. The air flow rate to 
either seeder was measured using a calibrated Brooks rotometer (Tube 
size R-6-15-B). The inlet pressure of the seeder was adjusted using a 
pressure regulator (Conoflow Model H40-XT-HXA). 

Particles were introduced into the main air stream using a 
vibrating, variable-speed, screw feeder (Vibra-Screw live bin feeder) 
with a 9.53 mm diameter screw. The screw discharged the particles 
into a tapered tube where they were allowed to drop downward into the 
main flow. The static pressures inside the particle feeder and the 
main air line were equalized at the particle entry position. A 45 m 
long, 3.18 mm inside diameter, tube was used to damp out any pressure 
oscillation caused by the introduction of the particles. Particle 
flowrates were measured by calibration of the feeder screw speed. 

The surroundings of the jet were also seeded with nominal 1 pm 
diameter aluminum oxide particles using a separate air supply line and 
a fluidized bed aerosol generator (TSI 9310). 

2.2.3 Swirl Generator 

A sketch of the swirl generator is illustrated in figure 2.4. 
Swirl was generated by introducing air tangentially into the injection 






tube through four identical 9.5 mm long slots located 90° apart. This 
approach for generating swirl has been previously reported by a number 
of investigators (6-11) and was found to yield reasonably good results 
during the present study. The slots were fed from a 62 mm inside 
diameter plenum, to insure equal flowrates through each slot. The 
leading edge of the slots was located 25 diameters upstream of the 
injection tube exit. The swirl number of the jet is varied by 
changing the flowrate of the tangential air flow. 

2.3 Instrumentation 

2.3.1 Single-Phase Jet Velocities 

The mean and fluctuating velocities of the single-phase jets were 
measured using a two-channel laser velocimeter (LV). A sketch of the 
arrangement appears in figure 2.5. Major components of the system are 
summarized in table 2.1. 

The output beam from a 2 W Argon-Ion laser was separated into 
green (514.5 nm) and blue (488 nm) beams. Each beam was then split 
Into two beams to provide a four-beam backward-scatter LV. This 
allowed simultaneous measurement of velocity and Reynolds stress in 
two directions. One beam from each color was frequency shifted (40 
MHz with electronic downshifting) to optimize the frequency range and 
to reduce fringe bias. Beam expansion was used to increase the 
signal-to-noi se ratio (SN) and to improve spatial resolution. The 
sendi ng/col lecti ng lens had a focal length of 762 mm. Theoretical 
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Figure 2.5. - Two-channel LV system. 
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parameters for this system, fringe spacing, measuring volume size, 
etc., are summarized in table 2.2. 

Detector signals were processed using two burst counters. An 
oscilloscope was used to monitor the detector and burst counter output 
signals. The counters were operated in the total burst mode, with a 
coincidence check, to insure that the signals obtained from both probe 
volumes were from the same seeding particle. Valid data rates were 
generally in the range of several hundred per second. The digital 
output of the counters was sent to a Mine 23 microcomputer, using a 
direct memory access board and a TSI interface (model 1998D). Data 
from the counters was processed using a Thermo-Systems computer 
program. The computer program included a residence time weighting 
scheme to minimize effects of velocity bias on the results. Results 
presented are based on at least 2400 velocity samples. Possible 
sources of error and experimental uncertainties for measurements of 
gas velocities in the single-phase jets are discussed in Appendix A. 
Experimental uncertainties (95 percent confidence) are as follows: 
mean streamwise velocities, less than 5 percent; mean angular 
velocities, typically less than 20 percent; velocity fluctuations, 
less than 5 percent; k, less than 12 percent; and the Reynolds stress, 
14 percent at the maximum Reynolds stress and proportionately higher 
elsewhere . 
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Table 2.2 - LV Theoretical Parameters 


Parameter 

Green 

Blue 


(514.5 nm) 

(488 nm) 

Single-phase Jets ( backseat ter ) : 



Fringe spacing, pirn 

4.76 

4.51 

Probe volume diameter, ptm 

102 

99 

Receiver measuring volume diameter, pun 

762 

762 

Probe volume length, pun 

1891 

1838 

Receiver measuring volume length, pim a 

103 

101 

Number of fringes 

22 

22 

Dispersed phase 30° off-axis (forward scatter); 



Fringe spacing, pun 

4.76 

4.51 

Probe volume diameter, yim 

102 

99 

Receiver measuring volume diameter, \jjvt 

685 

685 

Probe volume length, pun 

1891 

1838 

Receiver measuring volume length, pun® 

283 

277 

Number of fringes 

22 

22 


Calculated using depth of field limit. 

^Includes diameter of largest particle expected. 
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2 . 3.2 Particle - La den Jet V elocities 

2 . 3 . 2.1 Particle Velocities . - Mean and fluctuating velocities 
of the dispersed phase were measured using both a two-channel LV and a 
phase/Doppler anemometer. The LV system yielded particle velocities 
averaged over all particle sizes present while the phase/Doppler 
particle anemometer measured both particle size and velocity. 

The arrangement of the LV system used to measure particle 
velocities was similar to the LV system used to measure velocities in 
the single-phase jets, except that the receiving optics were moved to 
30° off-axis in the forward scatter direction. A receiving lens with 
a focal length of 602.4 mm and a clear aperture of 60 mm was used to 
collect scattered light from the particles. Estimated dimensions of 
the probe volume for this arrangement are also shown in table 2.2. To 
account for particles grazing the probe volume, the dimensions shown 
in table 2.2 were obtained by adding the diameter of the largest 
particle expected to the calculated theoretical diameter. Since 
actual probe volume size is a strong function of scattering particle 
size, laser power and photodetector gain setting, these values should 

not be strictly regarded as quantitative. 

The counters were operated in the total burst mode, however, a 
coincidence check was not required for this measurement (cross 
correlations were not required). The values obtained were number 
averages over all sizes of particles present in the flow. A residence 
time weighting correction was not applied or required, since only 
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number averages are meaningful for the dispersed phase. The 
continuous phase was not seeded with aluminum oxide particles for 
these measurements and laser power" and detector gain were also reduced 
to further insure that only signals from particles were processed. 

Valid data rates for particles varied considerably depending on where 
the flowfield was sampled. Each data point presented in the following 

represents at least 5000 samples. 

Possible errors and uncertainties for the particle velocity 
measurements are discussed in Appendix A. Experimental uncertainties 
(95 percent confidence) are as follows: u p , less than 5 percent, 

V , typically less than 10 percent (at x/d = 0.5); w p , less than 6 
percent at the maximum and proportionately higher elsewhere; and 
fluctuating particle velocities, less than 5 percent. 

The arrangement of the optical system for the phase/Doppler 
particle anemometer is illustrated in figure 2.6. Major components 
are summarized in table 2.3. The beam from a 2 W Argon-Ion laser was 
separated into green (514.5 nm) and blue (488 nm) beams. Only the 
green beam was used for these measurements. The beam was expanded by 
a factor of three before being split into two beams which were then 
focused and crossed with a 602.4 mm focal length lens. Scattered 
light was collected 30° off axis in the forward direction with a 495 mm 
focal length lens. The fringe spacing for this optical configuration 
was 6.3 pm. Estimated dimensions of the measuring volume are as 
diameter of 100 pm, and a length of 202 pm. These 
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Figure 2.6. - Single-channel phase/Ooppler particle anemometer. 



Table 2.3 - Summary of Components for Phase/Doppler Particle 

Anemometer 


Component 

Manufacturer 


Speci fi cation/function 

Argon-ion laser 

lexel 

95-2 

2.0 W, multiline mode; 1.3 mm 
(e~^ points) beam diameter 

Beam col 1 imator 

TSI 

9178 


Color separator 

TSI 

9105 


Beam expander 

Melles Griot 

09 LBM 001 

3.0 expansion ratio 

Polarization rotator 

TSI 

9102-12 


Beamspl i tter 

TSI 

9115-1 

50 mm beam spacing 

Transmitting lens 

TSI 

9119 

602.4 mm focal length 
2.34° half angle 

Receiving lens 

Aeromet ri cs 


495 mm focal length 

Receiver assembly 

Aerometri cs 

2100 

50 pm slit width 

240 mm focal length detector 

Processor 

Aerometri cs 

POP 3100 


Data acquistion/ 
processing 

IBM 

wm 
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measurements do not include the diameter of the largest particle. The 
collected light was focused onto a slit (50 pm wide by 1 mm long) 
using a 240 mm focal length lens. Three detectors, separated by fixed 
spacings, were used to detect the Doppler burst signals. The phase 
shift in the signal detected by the three detectors Is proportional to 
particle size and this shift can be used to determine particle size. 
Since this method also yields particle velocity in the same manner as 
a laser velocimeter, a simultaneous particle size and velocity 
measurement is obtained. Details of the theory and experiments to 
validate the instrument can be found in (76, 77). 

Data from the PDP3100 processor was processed using an IBM PC 
computer. The processor considered the entire Doppler burst for the 
data processing. The velocity data reported are number averages for 
each size group and no corrections were applied to account for velocity 
biasing. The actual number of samples considered for each velocity 
measurement varied depending upon position in the flowfield and 
particle size. Errors and uncertainties of the phase/Doppler 
measurements are discussed in Appendix A. The uncertainties (95 
percent confidence) are estimated to be less than 5 percent for mean 
and fluctuating streamwise velocities, similar to the single-phase jet 

measurements . 

232.2 Continuous Phase Velocities . - Mean and fluctuating 
velocities of the continuous phase in the presence of particles were 
measured using the single-channel phase/Doppler particle anemometer. 
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As previously discussed, this instrument measures particle size and 
velocity simultaneously. To measure continuous-phase velocities in 
the presence of particles, the continuous phase was seeded with 
nominal 1 pm diameter aluminum oxide particles and the corresponding 
velocity measurement for that size was assumed to represent the 
continuous phase. A reverse cyclone seeder, similar to that described 
by Glass and Kennedy (75), was used to introduce the aluminum oxide 
powder, since it provided Increased seeding levels over the smaller 
f luldized-bed seeder used previously for velocity measurements In the 
single-phase flows. 

By appropriate seeding of the continuous phase with the nominal 1 
pm diameter aluminum oxide particles, the phase/Doppler particle 
anemometer was used to measure continuous-phase velocities in the 
presence of the solid particles. An example of a typical data point 
taken with the instrument is shown in figure 2.7. In figure 2.7, two 
distinct size distributions are clearly evident. The smallest 
particles (~1.9 pm) are due to the continuous-phase seeding, while the 
larger particles correspond to the dispersed phase. The actual number 
of samples considered for each velocity measurement varied depending 
upon position in the flowfield, but at least 1000 samples were used 
for each continuous-phase data point. The velocity data reported are 
number averages and no corrections were applied to account for 
veloci ty biasing. 
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at a point 


Typical simultaneous size and velocity data taken 
in particle-laden flow. 


2.3.3 Particle Mass Flux 

The mass flux of the dispersed phase was measured using an 
isokinetic sampling probe. The probe used was similar to that 
reported by Szekely (78). A sketch of the probe appears in figure 
2.8. The probe was constructed using a Gelman 2220 stainless steel 
in-line filter holder. The probe is shown with a 2 mm inside diameter 
tip. Probe tips were interchangable and a tip with an inside diameter 
of 5 mm was also used. For the larger diameter probe, the height of 
the probe was reduced from 12.6 to 6.9 mm to increase the inlet 
diameter. Two probe inlet diameters were required in order to insure 
adequate spatial resolution and reasonable sampling times in all 
regions of the flow. The probe was rigidly mounted on a metal bar 
attached to a fixed optical table. The cage and injector tube 
assembly was traversed to sample various locations in the flow. 

A sketch of the particle mass-flux measurement system is shown in 
figure 2.9. A W.M. Welch Duo-Seal vacuum pump was used to withdraw 
samples i soki neti cal ly through the probe. Particles were collected on 
a Gelman membrane filter (No. 64679, 0.8 pm pore size) for a timed 
interval and weighed on a digital scale (Mettler PC 2000), accurate to 
0.01 g. The sampling interval was timed using an Adanac stopwatch. 
Sampling probe flow rates were controlled using a needle valve in the 
probe sample line and a second needle valve in a bypass line. 

Flowrates were measured using one of three calibrated rotometers (Cole 
Parmer FM102-05, FM082-03, FM034-39), arranged in parallel. Using 
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Figure 2.9. - Particle mass flux measurement flow system 



this sampling system, measured particle mass fluxes integrated across 
the jet were within ±10 percent of the calibrated feeder flow rate at 
all axial locations. This system was relatively insensitive to the 
effect of gas-sampling rates. Small variations in sampling rates 
above and below the mean isokinetic velocity had a negligible effect 
on the results. Experimental uncertainties for the particle mass flux 
measurements (95 percent confidence) are estimated to be less than 8 
percent, as discussed in Appendix A. 

2 . 4 Test Condi tions 

Test conditions are summarized in table 2 . 4 . Three single-phase 
and three particle-laden jets were studied. Swirl numbers were 
limited to 0.33 or less in order to avoid recirculation zones; 
therefore, the swirling flows could be classified as weakly swirling 
flows. A single loading ratio, defined as the ratio between the 
injected particle mass flow rate and the air mass flow rate, of 0.2 
was studied. Based on an average particle size of 30 pm, the spacing 
between particles at the exit of the tube was about 20 particle 
diameters. Particles with an SMD of 39 pm were used for all 
particle-laden flows. The standard deviation of the particle 
distribution was 15 pm. The average diameter of the particles was 30 
pm. The flows were fully turbulent, with initial Reynolds numbers of 
approximately 20,000 - based on the total jet momentum per unit mass 
flow rate and the tube diameter. 


39 




Structure measurements were made for x/d <30 for flows without 
swirl and x/d <20 for the swirling flows. For the particle-laden 
flows without swirl, radial profiles of flow properties were measured 
at x/d = 0.5, 5, 15, and 30. For the particle-laden swirling flows, 
radial profiles were measured at x/d = 0.5, 2, 5, 10, and 20. 
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CHAPTER III 


THEORETICAL CONSIDER ATIO NS 


3.1 General Description 

The analysis is limited to a steady, axi symmetr i c , dilute, 
solid-particle-laden, weakly swirling, turbulent jet in an infinite, 
stagnant media. The swirl number, calculated from equation (1.1), is 
restricted to values of approximately 0.5 or less so that adverse 
pressure gradients caused by the decay of angular velocity are not 
strong enough to induce a reversal of the axial velocity. The 
boundary-layer approximations are assumed to be valid; however, the 
radial pressure gradient, which is usually neglected in the 
boundary-layer analysis, is considered. The k-e turbulence model is 
used to provide closure since this approach has modest computational 
requirements. Effects of streamline curvature on the k-e turbulence 
model are considered. The injector exit Mach number is less than 0.3; 
therefore, kinetic energy and viscous dissipation of the mean flow are 
neglected with little error. 

Three methods of treating multiphase flow, typical of current 
practice, are considered, as follows: (1) locally homogeneous flow 

(LHF), where interphase transport rates are assumed to be infinitely 
fast and the flow can be treated like a single-phase, variable-density 


fluid; (2) deterministic separated flow (DSF), where finite 
interphase transport rates are considered, but the dispersed phase is 
assumed to interact only with the mean properties of the continuous 
phase e.g., particle/turbulence interactions are ignored; and (3) 
stochastic separated flow (SSF), where both finite interphase 
transport rates and effects of particle/turbulence interactions are 
considered, using random-walk methods. All three methods will be 
discussed in more detail in the following sections. 

3.2 Locally Homogeneous Flow 

3.2.1 Governing Equations 

3.2. 1.1 Baseline Version . - Since both phases move at the same 
velocity under the LHF approximation, this approach treats the flow as 

a single-phase turbulent fluid with density variations caused by changes 
in the concentration of the dispersed phase, even though the densities 
of each phase are constant. Following Bilger (79), Favre (mass)-averaged 
quantities are used rather than Reynolds (time)-averaged quantities. 

As discussed by Bilger (79) and Bradshaw, et al. (80), Favre averaging 
offers advantages for variable density flows because numerous terms 
involving density fluctuations are eliminated from the governing 
equations . 

Mean (mass averaged) quantities are found by solving governing 
equations for conservation of mass, momentum, and mixture fraction. 

The mixture fraction, f, was defined as the mass fraction of particles 
in the flow. If the diffusivities of both phases are the same, the 
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mixture fraction becomes a conserved scalar', as described in (79), and 


all scalar properties are only functions of mixture fraction. The 
equations are closed using the well-known, two-equation, k-c 
turbulence model . 

The general form of the governing equations in cylindrical 
coordinates is as follows: 


where 



<J> denotes a Favre-averaged variable, defined as follows: 


(3.1) 


P 

The source terms, S,, appearing in equation (3.1) are 

4 > 

summarized in table 3.1 along with empirical constants established by 
Jeng (81) for a variety of constant and variable density single-phase 
jets. The approach, and the empirical constants, however, are not 
very different from early proposals at Imperial College (74). The 
governing equation for radial momentum does not fit the form of 
equation (3.1) and is given later. 


Table 3.1. - Source Terms and Empirical Constants 
used i n Equation <3. 1 ) 



0.09 1.44 1.87 1.0 1.3 0.7 
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Equation (3.1) is classified as parabolic and was obtained using 
the usual boundary-layer assumptions. Axial gradients are assumed to 
be negligible compared to radial gradients, and turbulent correlations 
are all of the same order of magnitude. Thus, the radial dimensions 
of the flow are assumed to be much smaller than the axial dimensions, 
r / x < < 1 . Swirling flows exhibit fairly high spreading rates; 
therefore, the width of a swirling jet can be as large as 30 or 40 
percent of the axial dimension (11). Nevertheless, the boundary-layer 
equations have been applied to weakly swirling flows (swirl numbers up 
to about 0.5) (8, 10, 11, 16, 18, 22), and this approach is adopted 
here as well. 

In addition to the boundary-layer assumptions, a number of other 
assumptions are made in the derivation of equation (3.1). The flow is 
considered to be fully turbulent. Terms involving molecular viscous 
stress tensor fluctuations have been neglected, however, mean 
molecular stress is included even though it is considerably smaller 
than Also, a term involving the fluctuating density, velocity, 
and mean pressure gradient has been neglected in both the k and e 
equations in order to reduce the number of empirical constants. 

A consequence of the presence of angular velocity is that even 
though the standard boundary-layer assumptions are made, the radial 
momentum equation is still present as 


_~2 

9]3 _ pw 
8r r 


(3.3) 
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and cross-stream pressure gradients are not negligible. Also, because 
of the decay of angular velocity with axial distance, ttie source term 
3p/3x is included in the governing equation for axial momentum. 

Under the conserved scalar approach, the instantaneous particle 
concentration and density are only functions of mixture fraction. The 
mixture fraction is the mass fraction of particles in the flow, as 
fol lows : 


f = F -~- 1 V (3.4) 

a p 

Using the adiabatic mixing approximation from Shearer and Faeth (55), 
the density of the mixture can be found from the following equation: 

. ( p a - Op) f * -p 

L p a p P 

The Favre averaged mean value of a scalar can be found from the 
following equation (79,81): 


-1 


(3.5) 


$(f) 



4>(f )P(f)df 


(3.6) 


where 4>(f) is the state relationship for the property, i.e., 
equation (3.5) in this case. Time-averaged density can be found from 
from the following equation (81): 


1 

P 



- P(f )df 
p 


(3.7) 
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Since p~' is a linear function of f in the domain 0 < f < 1 , mean 
values are independent of the form of the probability density function 
(PDF), P(f), chosen and it is easily shown that p = p(f). Therefore, 
p . and f* can be substituted in equation (3.5) to give p . 
as a function of f. 

The turbulent viscosity appearing in equation (3.1) was 
calculated from k and e as follows (74): 



The source terms given in table 3.1 for k and e do not 
include effects of streamline curvature. This is considered in the 
next section. 

3. 3. 1.2 Streamline Curvature Version . - As discussed in Refs. 
82-88, streamline curvature has been shown to produce very large 
changes in the turbulence structure of shear layers. Bradshaw (82) 
presents a review of the effects of streamline curvature and discusses 
a basis for classifying shear flows as simple or complex. Shear 
flows, where there is only one significant rate-of-strai n component, 
are classified as simple while those subjected to extra rates of 
strain due to additional velocity gradients are classified as complex. 

Numerous modifications of the k-e turbulence model have been 
proposed to account for effects of streamline curvature. All of the 
modifications are based on the hypothesis that the destabilizing 
effect of swirl can be modeled through an increase in the length scale 
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of the turbulent eddies. Sharma (83) and Leschziner and Rodi (84) 
report calculations of curved flow using the k-c model, where one of 
the constants in the dissipation equation is made a function of the 
gradient Richardson number. Leschziner and Rodi (84) also reported 
results of calculations using a modified dissipation equation based 
upon a flux Richardson number. Hah and Lakshmi narayana (85) report 
predictions of turbulent wakes using turbulence closure models that 
were modified for effects of streamline curvature. A modification to 
the dissipation equation, obtained by simplification of an algebraic 
Reynolds-stress model, was reported by Leschziner and Rodi (86) and 
found to yield reasonably good results. Komori and Ueda (87) and 
Leschziner and Rodi (84), however, reported reasonably good results 
for calculations of strongly swirling free jets with standard k-e 
turbulence models as well. 

A recent review paper by Lakshmi narayana (88) discusses various 
turbulence models ranging from algebraic eddy viscosity models to 
models based on Reynolds-stress transport equations for complex 
flows. He states that "the field of turbulence modeling for complex 
flows is confusing and conflicting and that intuition and ad-hoc 
assumptions dominate the art of turbulence modeling in complex 
flows." Lakshmi narayana (88) concludes that the standard k-c model 
is not adequate for prediction of complex flows but modified k-e 
models are probably adequate for very mildly complex flows. Since 
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weakly swirling free jets are mildly complex, modifications to the 
k-c model are considered in the present study. 

Two modifications to the dissipation equations were evaluated 
during the present study. The first, discussed by Leschziner and Rodi 
(84), involves replacing C c2 in table 3.1 with the following 


expression 


C , = C ,<1 - 0.2 R.) 
z2 cZ 1 


where R. is a gradient Richardson number, which is defined as 
follows, for the present flows: 


k^ w_ 8(rw) 

i " 2 2 3r 

c r 


This curvature modification performed poorly during the present 
study and no predictions using it are presented in the following. 

The second modification of the dissipation equation also was 


taken from (84) and involves replacing one of the constants in table 
3.1 by a functional relationship. In this case, C el is replaced by 
the following expression: 

C*, = C ,(1 + 0.9 R.) (3.11) 

cl cl f 

where R^ is a flux Richardson number, defined in this case, as 
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3.2.2 Numerical Solution 


Equations (3.1) and (3.3) were solved using a modified version of 
the GENMIX computer code (89). A brief description of the formulation 
is given in Appendix B. Thirty-three cross-stream grid modes were 
used during all the calculations. For the nonswirling flows, 
streamwise step sizes were limited to 5 percent of the current flow 
width or an entrainment increase of 5 percent— whi chever was smaller. 
For the swirling flows, streamwise step sizes were limited to 2 
percent of the current flow width or an entrainment increase of 
2 percent - whichever was smaller. Computations performed with 
streamwise step size decreased by a factor of 2, and with 66 
cross-stream grid nodes, showed changes of less than 3 percent in the 
computed results at x/d » 20. 

Because of the formulation of the GENMIX algorithm, consideration 
of angular velocity required modifications of the solution procedure. 
Downstream values of r, which are required to solve equation (3.2) 
for the radial pressure distribution (in order to calculate the axial 
pressure gradient), can only be calculated after downstream values of 
u are known. Relatively accurate values of 3p/3x are required for 
the solution of the axial momentum equation. In order to deal with 
this problem, equation (3.1) for rw was solved first, to obtain 
downstream values of rw. Then equation (3.3) was integrated radially 
across the flowfield, to obtain radial pressures, using downstream 
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values of rw and upstream values of r. Then 9p/3x was calculated 
and then corrected based upon conservation of axial momentum flux. 
Using the corrected source term, 9p/3x, equation (3.1) was then solved 
for if, 7 , k, and e. This procedure has been used previously (10, 

11, 18). It preserves the marching character of the GENMIX solution 
procedure for the solution of the parabolic equations. For the 
present calculations, this procedure was found to conserve axial 
momentum within 2 percent. 

3.3 Separated Flow 

3.3.1 Continuous Phase 

The treatment of the continuous phase in the separated-f low 

analysis is similar to the LHF model, except that additional source 

terms, S due to interphase transport are included in the 
P9 

governing equations. The flows considered here are very dilute; 
therefore, effects of the dispersed phase on turbulence quantities are 
ignored for the present. In addition, since isothermal, solid, 
particle-laden flows are considered, the only source term involves the 
exchange of momentum between the continuous and the dispersed phase. 
Because the density of the continuous phase is constant and its volume 
fraction is nearly unity, a solution for f is no longer required. 
Finally, in this case, Favre and Reynolds averaged quantities are 
identical . 

Both separated-f low analyses involve dividing the dispersed phase 
into n groups and tracking the trajectory of each group through the 
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flow field. In order to obtain the momentum exchange source term, the 
continuous-phase flow field is divided into computational cells and 
the net change in momentum of each particle group as it enters and 
leaves a computational cell is computed. This is referred to as the 
"Particle Source in Cell Approach" (PSIC) described by Sharma and 
Crowe (59). The exact forms of the source terms are given in table 
3.2. In order to preserve the marching character of the calculation, 
no source term is calculated for the radial direction. This can be 
neglected with little error since the radial momentum equation is only 
used to calculate the axial pressure gradient which is later corrected. 

Particle source terms could also be included in table 3.2 for 
both k and e, to represent the effects of the particles on 
turbulence properties. As discussed by A1 Taweel and Landau (90), 
turbulence intensities of the continuous phase can be reduced as a 
result of the presence of particles. They concluded that the 
magnitude of this damping, termed turbulence modulation, increases 
with increasing particle loading and decreases with increasing 
particle diameter. There have been a number of attempts to account 
for the influence of the dispersed phase on turbulence properties 
using the k-e turbulence model. Mostafa and Elghobashi (91) derived 
a two-equation turbulence model for two-phase flows that attempts to 
account for the additional energy dissipation due to the presence of 
particles. Mostafa and Mongia (92) evaluated a simplified version of 
this model which contained one additional constant and found that it 
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Table 3.2. - Particle Source Terms in 
the Separated Flow Analysis 


* 


s 

P<J> 

1 


0 

u 

n 

E 

J i = l 

n m [u \ - (u . ) 

i P PV • \ P 1 / f 

i n out 

' J 

rw 

-1 A 

v j 2 

J i = l 

f — 

n m /rw . \ - /rw 

1 P 1 pi 'in ' p, 'out 
w J J 

k 


0 

e 


0 


54 



yielded slightly better predictions than the standard k-c model when 
the constant was optimized for that particular flow considered. 

Shuen, et al . (40), Shuen (41), and Zhang, et al . (42) also presented 
a model for turbulence modulation based on the k-c model. 

Comparisons with measurements were generally better than the 
single-phase k-c model, however, an additional constant was also 
introduced which could not be easily evaluated. Since the loading 
ratio of particles is relatively low ( 0 . 2 ) in the present study, 
turbulence modulation was not considered in the theory. 

3.3.2 Dispersed Phase 

The di spersed-phase properties are obtained by solving the 
Lagrangian equation of motion for the particles, assuming that the 
particles can be approximated as spheres. The general form of the 
equation (the B-B-0 equation that includes effects studied by Basset, 
Boussinesq, and Oseen), after neglecting effects of particle rotation, 
can be written as follows (27): 


it . 3 u d ir , 2 fs i *♦ I /"* x it , 3 3p 
7 d p . 4 / = o d pC n u - u (u - u ) - t d 
6 p K p dt 8 p M D 1 p 1 p 6 p 


tr , 3 d + \ 3 . 2 , x 1 12 

* \2 d p p 5 T <u ‘ V * 2 d p <w> 


(d/d?(u - 3 >) 

Jt w-T « * F s 

- ,t po (t - V 


(3.13) 

where the time derivative is taken following the motion of the particle 

hT~ = If + ~ (3.14) 

dtp 3t p 


55 


The term on the left-hand side of equation (3.13) represents the 
inertia of the particle sphere. Taken in order, the terms on the 
right-hand side of equation (3.13) represent: the drag force on the 

sphere, which conventionally includes both skin friction and form 
drag; the force on the sphere due to static pressure gradients in the 
flow; the force on the sphere due to the inertia of fluid displaced by 
its motion, which is often called the virtual mass term; the Basset 
term, which allows for effects of the deviation of the flow from a 
steady flow pattern around the sphere; and the external or body force 
term, e.g., the force due to gravity. 

Assumptions made to calculate particle trajectories were as 
follows: dilute particle-laden flow with drag equivalent to a single 

particle in an unbounded environment; particle collisions neglected; 
drag treated empirically, assuming quasisteady flow for spherical 
particles; and since Pp /p>200, effects of static pressure 
gradients, virtual mass, Basset forces, Magnus forces, etc., can be 
neglected with little error. The remaining assumptions are typical of 
separated-f low models of dilute particle-laden flows and are described 
more completely in Refs. 54 and 57. 

After adopting these assumptions, equation (3.13) can be greatly 
simplified with the result that the position and velocity of each 
particle group can be found by integrating: 

dx . 

= u ., i = 1,2,3 (3.15) 

dt pi 
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The drag coefficient was calculated from a standard empirical 
correlation for solid spheres, recommended by Faeth (62-54), as 
fol lows : 


C 


D 


_ 24_ 
= 



p £ '000 


C Q = 0.44, Re p > 1000 


(3.17) 


Equation (3.16) was integrated using a second-order algorithm similar 
to Shuen (41 ) . 

3.3.3 Deterministic Separated Flow 

For the deterministic separated flow approach, the dispersed 
phase is assumed to interact only with mean properties of the 
continuous phase. Since mean properties are only considered for the 
deterministic separated flow model, all terms in equation (3.16) are 
considered to be time averaged and the local mean velocity is used as 
the ambient velocity for each particle group. 

The coupled parabolic equations, (3.1), (3.3), and (3.16) were 
solved in an iterative fashion. Equations (3.1) and (3.3) were solved 
first with the particle source term from the previous step. Then 
equation (3.16) was solved to determine particle trajectories and to 
update source terms for the continuous phase. At least 1400 particle 
groups were tracked during these calculations in order to find 
statistically significant flow properties. 
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3.3. 4 Stochastic Separated Flow 

The stochastic separated flow approach involves finding the 
motion of a statistically significant sample of particles as they 
leave the injector and encounter a succession of turbulent eddies. 
Treatment of the continuous phase is identical to the deterministic 
separated flow model. Equations (3.15) and (3.16) are also solved for 
particle motion, however, instantaneous properties are used rather 
than mean continuous-phase properties. The method of specifying 
instantaneous eddy properties and treating particle/eddy interactions 
follows a proposal of Gosman and Ioannides (66), as modified by Shuen 
(41 ) . 

Important properties of the stochastic separated flow computation 
are the physical properties of each eddy and the interaction time 
between a particle group and particular eddy. The velocity within 
each eddy is assumed to be constant but eddy properties are assumed to 
vary randomly from one eddy to the next. Eddy velocities are found 
using Monte Carlo methods, after defining a probability density 
function (PDF) for each component of velocity. Velocity fluctuations 
are assumed to be isotropic, with a Gaussian PDF having a standard 
deviation of (2 k/3) and mean values u, v, and w. The cumulative 
distribution function (CDF) of each velocity component is constructed 
and randomly sampled. This involves choosing three random numbers in 
the range 0-1 with a random number generator and then finding the 
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three velocities at these values of the CDF. This procedure assures 
random selection of velocities in a manner that satisfies the PDF of 
velocity. 

Particle groups are assumed to interact with a particular eddy fot 
a time either as long as the eddy lifetime or the time required for a 
particle to traverse the eddy, whichever is shorter. Following Shuen 
(41), characteristic eddy sizes and lifetimes are specified as follows: 



Particles and eddies are assumed to interact as long as the time 

of interaction and the relative displacement of the particle and the 

eddy (from the start of the interaction) are both less than t g and 

L When the interaction ends by the L criterion, the particles 
e ' e 

have traversed the eddy. Ending the interaction by the t^ criterion 
implies that the eddy has captured the particle. 

The remaining computations are similar to the deterministic 
separated flow model . The random-walk calculations for each particle 
group, however, required a larger number of particle groups to obtain 
statistically significant results: at least 7000 particle groups were 

employed during the present computations. 
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CHAPTER IV 


RESULTS AND DISCUSSION 


4.1 Single-Phase Jets 

Single-phase jets were studied initially in order to assess the 
suitability of the experimental configuration, to provide baseline 
results for comparison with measurements for the particle-laden jets, 
and to assess the capability of the continuous-phase model without the 
complications resulting from the presence of a dispersed phase. 
Experimental results obtained during the present study are discussed 
first. Predictions are then compared with measurements obtained both 
during the present study, as well as other measurements in swirling 
jets taken from the literature. All data obtained during the present 
study are tabulated in Appendix C. 

4.1.1 Experimental Results 

Measured properties of the present single-phase jets are 
illustrated in figures 4.1 to 4.13. Based on the measurements 
illustrated in figures 4.1, 4.2, and 4.3, for S = 0, 0.19, and 0.33, 
respectively, it is clear that swirl has a large effect on streamwise 
properties along the flow axis. The rate of decay of the axial 
velocity with streamwise distance, seen in figure 4.1, increases as 
swirl number increases. For example, the mean centerline axial 
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Figure 4.2. - Axial variation of centerline turbulence kinetic 
energy with streamwise distance for the single-phase nonswirling 
and swirling jets (S = 0, 0.19, and 0.33). 





Figure 4.3. - Axial variation of maximum angular velocity with 
streamwise distance for the single-phase swirling jets 
(S = 0.19 and 0.33). 
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Figure 4.4. - Radial variation of flow properties at x/d = 0.5 for 
the single-phase nonswirling jet. 
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Figure 4.10. - Radial variation of flow properties at x/d = 20 for 
the single-phase swirling jet (S = 0.19). 
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velocity decays to one-half its original value (measured at x/d = 
0.5) in approximately 3.5 injector diameters for S = 0.33, while 
approximately 13 injector diameters are required for the jet without 
swirl. Turbulence kinetic energy near the injector exit increases 
dramatically with increasing swirl (see figure 4.2). At x/d = 0.5, 
normalized k c is approximately 25 times higher for 5 = 0.33 than 
for the nonswirling jet. For the swirling jets, the development of 
k c with streamwise distance is strongly influenced by the swirl 
number. For S = 0.19, shown in figure 4.2, normalized k c begins 
to increase near the injector exit and rises monotonical ly to an 
asymptotic value on the order of 0.10. For S = 0.33, also shown in 
figure 4.2, normalized k c increases even more sharply near the 
injector exit, reaching a peak of approximately 0.22 at x/d = 2, and 
then decays to approximately the same asymptotic value as S = 0.19. 
Angular velocity, see figure 4.3, decays quite rapidly with streamwise 
distance. For S = 0.33, angular velocity has decayed to one-half its 
original value (at x/d = 0.5) at approximately x/d = 2. Increasing 
swirl increases the rate of angular velocity decay. 

Radial profiles of measured properties in the single-phase jets 
are illustrated in figures 4.4 to 4.6 for S = 0, figures 4.7 to 4.10 
for S = 0.19, and figures 4.11 to 4.13 for S = 0.33. At x/d = 0.5, 
peak values of k are increased by a factor of two for S = 0.33, 
when compared with the no-swirl case as shown in figures 4.4 and 
4.11. Similar results are observed for Reynolds stress. The results 



illustrated in figures 4.4, 4.7, and 4.11 also indicate that, at 
x/d = 0.5, measured values of velocity fluctuations increase as the 
swirl number is increased, as expected. Even at x/d = 0.5, 
increasing the swirl number increased the jet width; this can be seen 
by comparing figure 4.4 for S = 0, and figure 4.11 for S = 0.33. 

Radial profiles of measured properties for the swirling jets at 
x/d = 5 are illustrated in figures 4.8 and 4.12. As expected, the jet 
width for the higher swirl number jet (S = 0.33) is larger than the 
S = 0.19 jet. Measured values of fluctuating quantities also increase 
with increasing swirl number. For S = 0.33, shown in figure 4.12, 
measured values of all components of the velocity fluctuations are 
approximately equal at x/d = 5. For S = 0.19, however, the 
streamwise velocity fluctuations are slightly higher than the other 
two components at this axial location. 

Radial profiles of measured properties for the single-phase jets 
with swirl at x/d = 10 appear in figures 4.9 and 4.13 for S = 0.19 
and S = 0.33, respectively. Angular velocities have decayed to 
approximately 0.4 m/s for both swirling jets at this location, and as 
a result, the structure of both jets is similar. However, the jet 
width is still wider for the higher swirl number jet than for the 
lower swirl number jet. 

Radial profiles of measured properties for the S = 0.19 swirling 
jet at x/d = 20 are illustrated in figure 4.10. The mean angular 
velocity has decayed to a negligible value and measured values are 
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similar to a jet without swirl at this location. For comparison, 
radial profiles of the nonswirling jet at x/d = 15 and 30 appear in 
figures 4.5 and 4.6, respectively. 


4.1.2 Single-Phase Jet Pred ict ions 

In order to initiate evaluation of predictions, calculations were 
completed for comparison with representative measurements for weakly 
swirling single-phase free jets, taken from the literature. Only 
measurements where initial conditions were known reasonably well were 
considered. Swirl numbers for these flows ranged from 0.25 to 0.5. 

Initial conditions for the predictions are taken at the 
measurement location nearest to the injector exit. Starting values 
of c were determined from the definition of a turbulent length 
scale, as follows: 


= c 

o p L 


(4.1) 


where L was chosen as a fraction of the initial jet half-width to 
provide reasonably good agreement with initial streamwise changes of 
k. The value used for L was on the order of 0.3 of the initial 
half-width of the flow. 

Results using the measurements reported by Morse (11) are 
illustrated in figures 4.14 to 4.17. Morse (11) obtained measurements 
using a single hot-wire probe for two weakly swirling free jets having 
swirl numbers of 0.25 and 0.35. Initial flow properties were measured 
at x/d =0.5. Predictions using the standard k - c turbulence 
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model and the modification for streamline curvature based on the flux 
Richardson number (equation (3.12)) are shown on the figures. 

The streamwise variation of flow properties for S = 0.25 is 
illustrated in figure 4.14. For these results, the approach using a 
curvature correction gave better agreement with the measurements than 
the standard k - e turbulence model, however, the differences 
between the two are not very large. The curvature modification, 
however, caused predicted values of k along the axis to be 
significantly overestimated. Measurements of k c reported by Morse 
(11), however, are somewhat lower than values obtained during the 
present study, which may be a factor in this apparent deficiency of 
the curvature correction approach. Radial profiles at x/d = 10 for 
measurements with S = 0.25 reported by Morse (11), are shown in 
figure 4.15. Again, the predictions are reasonably good. The 
standard k - c turbulence model somewhat underestimates the jet 
width at x/d = 10 based on the streamwise mean velocities while use 
of the flux Richardson correction causes the width of the jet to be 
slightly overestimated. The comparison between predicted and measured 
values of k and the Reynolds stress, illustrated in figure 4.15, is 
similar to the mean streamwise velocity. 

Measurements and predictions for a jet having S = 0.35, reported 
by Morse (11), are illustrated in figures 4.16 and 4.17. For the 
higher swirl number flow, predictions using the curvature correction 
are in better agreement with measurements of the decay of mean axial 



and maximum angular velocity. The development of k with streamwise 
distance is again overestimated using the curvature correction, 
however, measured values of k reported by Morse (II) were somewhat 
lower than those obtained during the present study. Normalized values 
for the decay of mean centerline axial and maximum angular velocity 
reported by Morse (11) show good agreement with those obtained during 
the present study. Radial profiles for the S = 0.35 jet measurements 
reported by Morse (11) at x/d = 10 are shown in figure 4.17. Again, 
predictions using the curvature correction show better agreement with 
measurements than the standard k - c model, although the jet width 
at x/d = 10 is slightly underestimated, even with the curvature 
correction. Morse ( 11 ) states, however, that edge values of angular 
velocity are not reliable and should be used with caution. The 
agreement between the radial profiles at x/d = 10 reported by Morse 
( 11 ), shown in figure 4 . 17 , and measurements obtained in the present 
study for similar conditions, figure 4 . 13 , is reasonably good. 

Measurements at a higher swirl number, S = 0.5, are reported by 
Sislian and Cusworth (5). These measurements were obtained using a 
single-channel LV system. Initial conditions were measured very 
close to the injector exit, at x/d = 0.125. Measurements were only 
made near the jet exit, ending at x/d = 5. The variation of 
streamwise flow properties is shown in figure 4.18 for this flow, 

Mean centerline axial and maximum angular velocities are predicted 
reasonably well for x/d > 3, using the curvature correction. 









Predicted rates of decay are overestimated, however, very close to the 
injector exit. Since the swirl number is relatively high for this 
flow, it may not be well represented by parabolic governing equations. 

Measurements and predictions of radial flow profiles at x/d = 4 
are shown in figure 4.19, for the measurements of Sislian and Cusworth 
(5). The agreement between the predictions and measurements is quite 
good for the approach using the curvature correction. No radial 
profiles at streamwise distances greater than x/d = 4 were reported. 

Comparisons between the measurements and predictions obtained 
during the present study are illustrated in figures 4.20 to 4.22 for 
the single-phase jets without swirl. The predictions generally are in 
reasonably good agreement with the measurements for the nonswirling 
jets. For example, predictions of the rate of decay of mean 
centerline axial velocities with streamwise distance and radial 
profiles of axial velocity, both show good agreement with measurements. 
However, at x/d = 15, illustrated in figure 4.21, predictions of k 
and Reynolds stresses are slightly higher than measurements, 
particularly near the flow axis. In contrast, the jet width at this 
position is predicted quite well. Farther downstream, predictions of 
k and Reynolds stress are in better agreement with measurements, see 
the findings at x/d = 30, illustrated in figure 4.22. 

Predictions and measurements are illustrated in figures 4.23 to 
4.26 for the S = 0.19 weakly swirling jet of the present study. As 
shown in figure 4.23 for flow properties along the axis, the approach 











Figure 4.20. - Predicted and measured variation of flow properties 
in the streamwise direction for the single-phase nonswirling jet 































using the streamline curvature correction based on the flux Richardson 
number gave better results than the standard k - c turbulence 
model. Predictions and measurements of radial flow properties at x/d 
= 5, 10, and 20 are shown in figures 4.24, 4.25, and 4.26, 
respectively. Again, predictions using the curvature correction are 
in better agreement with measurements than the standard k - c 
turbulence model. The curvature modification always increases levels 
of k across the entire width of the jet. This increases both the 
jet width and rate of decay of mean velocity components with 
streamwise distance. As expected, as the angular velocity decreases, 
differences between predictions for turbulence models with and without 
the curvature modification are decreased as well. 

Results are not shown for the curvature modification based on the 
gradient Richardson number (equation (3.10)). This modification 
produces a turbulence damping effect over regions of the jet where 
8rw/3r is positive. Levels of k are only increased in the edge 
regions of the jet where 3rw/3r is negative. This is clearly not 
correct for the swirling free jets considered during the present 
study. In contrast, the curvature modification based on the flux 
Richardson number tends to increase turbulence levels regardless of 
the sign of 3rwV3r. 

Comparisons between measurements and predictions for S = 0.33, 
from the present study, are illustrated in figures 4.27 to 4.29. Flow 
properties along the axis are illustrated in figure 4.27 and 
demonstrate that predictions are reasonably good, especially when 
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using the curvature correction. The sharp increase in k at x/d - 2 
is not predicted very well, however, even with the streamline 
curvature modification. The variation of mean centerline axial and 
maximum angular velocity with streamwise distance is well piedicted. 
Predicted and measured radial variations of flow properties for the 
S = 0.33 jet of the present study at x/d = 5 and 10 are shown in 
figures 4.28 and 4.29, respectively. Again, the curvature 
modification shows better agreement with measurements than the 
standard k - e turbulence model. Predictions are reasonably good at 
both axial locations. 

The predicted static pressure variation with streamwise distance 
for the S = 0.19 and 0.33 swirling single-phase jets is illustrated 
in figure 4.30. Static pressures along the centerline are slightly 
below atmospheric pressure. For the higher swirl number, S = 0.33, 
the maximum pressure difference is 70 Pa. The static pressure along 
the axis rapidly approaches atmospheric pressure since angular 
velocity decays rapidly. Predictions illustrated in figure 4.30 were 
obtained with the flux Richardson number modification to the k - c 
turbulence model. Predicted static pressures using the standard 
k - c model were similar. 

Based upon the single-phase results, the flow structure of the 
single-phase weakly swirling jets appears to be reasonably well 
predicted with turbulence models that were considered during this 
study. The modification of the dissipation equation to account for 
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effects of streamline curvature, based on the flux Richardson numbet , 
generally improved the predictions. 

4 . 2 Part i c 1 e -L a de n Jets 

4 . 2 . 1 Nc_-Swi r 1 Condi t ion s 

Measurements and predictions, using the three two-phase flow 
models for the nonswirling jet are discussed in this section. Initial 
conditions for the particle-laden jets were measured at x/d = 0.5: 
the results of these measurements are illustrated in figures 4.31 and 
4.32 for, S = 0. 

Figure 4.31 illustrates number-averaged particle velocities which 
are averaged over all size groups, while figure 4.32 illustrates mean 
particle axial velocities for particle diameters of 23, 43, and 63 
pm. From figure 4.32, it is evident that at x/d = 0.5, the larger 
particles are moving slower than the smaller particles at the center 
of the jet, but their radial profile is flatter than the smaller 
particles so the larger particles are moving faster near the edge of 
the jet. Differences in axial velocity between size groups are quite 
large near the edge of the jet. 

A single loading ratio of 0.2 was used for all test conditions 
(summarized in table 2.4) and the particles had an SMD of 39 pm with a 
standard deviation of 15 pm. Because of instrument limitations, only 
mean and fluctuating axial velocities of the gas phase could be 
measured, as discussed in the Experimental Methods section. 

Therefore, it was necessary to estimate initial values of turbulence 
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kinetic energy for the gas phase. For the predictions shown, initial 

values of k were assumed to be the same as the single-phase flows. 

2 1/2 

Measured values of (u 1 ) for the gas phase were approximately 
20 percent lower across the entire jet width for the particle-laden 
jets (at the initial condition of x/d = 0.5) than for the 
corresponding single-phase jets. Predictions showed that reductions 
of k of 20 percent caused negligible changes in flow properties, 
however, except very close to the injector. Initial values of c for 
the particle-laden jet predictions were also unchanged from the 
single-phase jets. For predictions using the separated-f low models, 
seven particle size groups ranging in size from 12.5 to 72.5 microns 
in 10 micron increments were tracked through the flowfield. Measured 
mean and fluctuating particle axial velocities at x/d = 0.5 from the 
phase/Doppler particle anemometer were used as initial conditions. 
Initial values of mean and fluctuating particle radia 1 velocities were 
taken from LV measurements at x/d = 0.5, which were number averaged 
over all particle sizes. Fluctuating gas-phase axial velocities are 
not shown in figure 4.32, however, they, along with all data taken 
during the present study, are tabulated in Appendix C. 

The predicted and measured variation of axial velocities in the 
streamwise direction for the nonswirling, particle-laden jet are 
illustrated in figure 4.33. Measured particle velocities reported in 
figure 4.33 are number averaged over all particle sizes. For 
comparison with these measurements, predictions from the 
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separated-f low models were number averaged over all seven size groups 
as well. By comparing figures 4.1 and 4.33, it is evident that axial 
velocity decays more slowly for the particle-laden flow than for the 
single-phase flow due to momentum exchange from the particles. 

Because the particle loading ratio is relatively low, predictions from 
both the locally homogeneous flow and stochastic separated flow models 
are nearly identical and also show reasonably good agreement with the 
experimental measurements. For the gas-phase axial velocity, 
deterministic separated flow model predictions were identical to those 
from the stochastic separated flow model and are not shown. 

Predictions using the locally homogeneous flow and separated-f low 
models for particle axial velocities show distinct differences. The 
neglect of particle inertia, illustrated by predictions from the 
locally homogeneous flow model, overestimates the rate of decay of 
particle axial velocity. Predictions from both the stochastic and 
deterministic models are similar and show reasonably good agreement 
with measurements for particle axial velocity decay. 

Radial profiles of flow properties, number averaged over all 
particle sizes, for the particle-laden jets without swirl, are 
illustrated in figures 4.34 to 4.36. Measurements and predictions are 
illustrated for x/d = 5, 15, and 30. Gas-phase jet widths are 
slightly overestimated for the particle-laden flows at x/d = 15 and 
30, see figures 4.35 and 4.36. The overestimation of gas-phase flow 
width is probably caused by turbulence modulation due to the 
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particles. Even though the loading ratio was relatively low, values 
of < u ' ) for the particle-laden jets are reduced compared to 
the single-phase jets. The reductions in fluctuating velocities 
reduce turbulent mixing rates and, thus, the jet width. Particle 
axial velocity profiles are reasonably well predicted by all the 
models, however, since the locally homogeneous flow model 
overestimates the rate of axial velocity decay with streamwise 
distance, the values of axial velocity predicted by this approach are 
generally lower than the measurements. 

Radial profiles of particle mass flux illustrate the different 
physical assumptions embodied in the three models. As illustrated in 
figures 4.35 and 4.36, the no-slip assumption of the locally 
homogeneous flow model causes the radial dispersion of the particles 
to be overestimated. At x/d < 15, both separated-f low models give 
similar predictions of particle mass flux. However, at larger 
streamwise distances, see figure 4.36, the neglect of turbulent 
dispersion of the particles causes the deterministic separated flow 
model to underestimate the spread of the particles. Since mean 
gas-phase radial velocities are quite small at these axial distances, 
turbulent dispersion is the only mechanism available for radial spread 
of the particles. Only the stochastic separated flow model, which 
accounts for both particle inertia and turbulent disperison of the 
particles, correctly predicts the radial distribution of particle mass 
flux at x/d = 30 (see figure 4.36). 


The predictions of fluctuating particle properties for the 

stochastic separated flow model, also illustrated in figures 4.34 to 

_t 2 1/2 

4.36, are reasonably good, however, (u ) is slightly 

underestimated at x/d = 30. This is probably caused by the 

assumption of isotropic velocity fluctuations, which causes streamwise 

continuous-phase velocity changes experienced by the particle to be 

~ r 2 1 / 2 

underestimated. In particular, measured values of (u ) were 

' 21/2 ' 21/2 

always greater than (v ) and (w ) for the single-phase 
jet and it is expected that this behavior should be similar for the 
particle-laden jet. 

These results for the particle-laden jets without swirl are 
similar to those previously reported by other investigators for 
similar flows (38 to 41 ) . 

Measurements and predictions of particle mean axial velocity for 
23, 43, and 63 pm particles at x/d = 5, 15, and 30 are illustrated in 
figure 4.37 for the nonswirling jet. Only predictions from the 
stochastic separated flow model are presented. As illustrated in 
figure 4.37, decay of axial velocity of the larger particles is slower 
than the smaller particles. By x/d = 15, the 63 micron particles 
have the highest velocity. Predictions from the stochastic separated 
flow model are in reasonably good agreement with measurements and 
correctly show the general trends. 

4.2.2 Swirling Conditions 

Measurements and predictions for the particle-laden swirling jets 
are discussed in this section. Initial conditions for the 
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particle-laden swirling jets were measured at x/d - 0.5. Jets having 
swirl numbers of S = 0.16 and 0.3 were studied. A single loading 
ratio of 0.2 was also used for the swirling flows. Test conditions 
are summarized in table 2.4. As previously discussed, IV particle 
velocity measurements were number averaged over all particle sizes. 

In addition, simultaneous particle sizes and axial velocities were 
measured using the phase/Doppler particle anemometer. 

Measurements of flow properties, averaged over all particle sizes, 
are illustrated in figures 4.38 and figure 4.39 for the particle-laden 
swirling jets. Mean axial velocities of 23, 43, and 63 pm particles 
at x/d = 0.5, are illustrated in figure 4.40. It is apparent that 
the presence of swirl has a large effect on the particle-laden jets, 
see figures 4.38 and 4.39. Compared to the jet without swirl, the 
maximum particle mass flux has shifted outward from the centerline of 
the jet to r/x = 0.8. Particle velocity fluctuations also increase 
with increasing swirl. In addition, increasing swirl also increases 
the variation of mean axial velocity with particle size, as 
illustrated in figure 4.40. At x/d = 0.5, larger particles are 
moving at a lower axial velocity than smaller particles throughout 
most of the jet. The radial velocity profile of the larger particles 
is much flatter than the smaller particles: as a result, the 

velocities of the larger particles are greater than the smaller 
particles near the edge of the jets. Fluctuating particle axial 
velocities for each size group are not shown in the figures but are 
tabulated in Appendix C. 
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Figure 4.39. - The number-averaged initial conditions for the 
particle-laden swirling jet (S = 0.3) — radial variation of 
flow properties at x/d = 0.5. 
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Figure 4.40. - Radial variation of mean axial velocity for 
23, 43, and 63 micron particles at x/d = 0.5 for the 
particle-laden swirling jets (S = 0.16 and 0.3). 


As discussed earlier for the particle-laden jets, only mean and 
fluctuating axial velocities of the gas phase could be measured. 

Thus, it was necessary to estimate initial values of the gas-phase 
angular velocity and the turbulence kinetic energy in order to 
initiate calculations. The initial angular velocity of the continuous 
phase was estimated by subtracting the measured initial particle-phase 
angular momentum from the values obtained for the single-phase flows. 
For both swirling flows, particle-phase angular momentum was 
approximately 10 percent of the measured single-phase angular 
momentum. Initial values of k and also e were assumed to be the 
same as the single-phase flows. Separated-f low predictions were made 
using seven particle size groups having diameters ranging from 12.5 to 
72.5 pm, in 10 pm increments. Initial values of particle mean and 
fluctuating axial velocities for each size group were obtained from 
phase/Doppler measurements at x/d = 0.5. Initial values of particle 
radial and angular velocities were assumed to be identical for all 
size groups, and were obtained from number-averaged LV measurements. 
Sensitivity of the predictions to initial values is discussed in a 
later section. 

Predicted and measured flow properties in the streamwise 
direction are illustrated in figures 4.41 to 4.43 for the swirling 
particle-laden jets. Only the standard k - e turbulence model was 
employed for predictions of the particle-laden flows reported in the 
following. Reasons for not considering a correction for streamline 
curvature are discussed later. 


125 











6 

.2 
0 

1.0 
.8 
6 
M 

.2 

0 

ure 4.43. - Stre 
laximum angular v 
S = 0.16 and 0.3 






Predicted and measured gas-phase axial velocities are shown in 
figure 4.41 for swirl numbers of 0.16 and 0.3. As illustrated in 
figure 4.41, increasing swirl increases the rate of decay of axial 
velocity, similar to single-phase flows. Predictions from the locally 
homogeneous flow and stochastic separated flow models are nearly 
identical and show reasonably good agreement with measurements for 
both swirling flows. Again, predictions of gas-phase axial velocity 
using the deterministic separated flow model are nearly identical to 
predictions using the stochastic separated flow model and are not 
shown in the figure. 

Particle mean axial and maximum angular velocities, number 
averaged over all particle sizes, are plotted as a function of 
streamwise distance in figures 4.42 and 4.43, respectively. As 
expected, neglecting slip between the phases causes predictions from 
the locally homogeneous flow model to overestimate the rate of decay 
of particle velocities. Differences between predictions from the 
deterministic and stochastic separated flow models were small: both 

show good agreement with measurements. At the initial condition (x/d 
= 0.5), axial velocities of the particles were lower than the 
continuous phase except near the edge of the jet for both swirling 
flows, see figures 4.38 and 4.39. Particle axial velocities initially 
increase, due to momentum exchange from the continuous phase, before 
beginning to decay. As can be seen from figure 4.42, both 
separated-f low models correctly predict this behavior. 
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Radial profiles of flow properties for the particle-laden 
swirling jets are illustrated in figures 4.44 to 4.53. Velocity 
measurements, number averaged over all particle sizes, are illustrated 
in figures 4.44 to 4.51. Particle mean axial velocities for three 
size groups are illustrated in figures 4.52 and 4.53 for the two 
swirling flows. 

Measurements and predictions of radial flow profiles at x/d = 2 
are illustrated in figures 4.44 and 4.45 for S = 0.16 and S = 0.3, 
respectively. Similar to the results observed for the nonswirling 
particle-laden flow, predictions of gas-phase velocities using the 
locally homogeneous flow and stochastic separated flow models are 
nearly identical and both slightly overestimate the gas-phase flow 
width. Again, the probable cause is turbulence modulation by the 
particles, which was not considered by the present models. All 
predictions reported here employed the standard k - c turbulence 
model. For single-phase jets with swirl, the jet width was always 
underestimated using the standard k - e turbulence model and 
streamline curvature modifications were introduced which increased the 
predicted jet width. Because of the turbulence modulation by the 
particles, predictions with the streamline curvature modification 
showed even poorer agreement with the measurements than the standard 
k - e model and are not illustrated in the figures. At x/d = 2, the 
stochastic separated flow model overestimates particle velocities for 
r/x greater than about 0.2. Predictions of particle axial velocity 




















.3 , — 



O o o 


cm 
__ c 
!=> 


• 1 r- 


o >- 


u 



cm 

Q. 

I > 


.2 

0 


u 



rsj 


.2 

.i« 


0 


Figure 4.48. - Radial var 
x/d = 10 for the parti c 



















from the deterministic separated flow model show the best agreement 
with the measurements at x/d = 2. The locally homogeneous flow model 
underestimates particle axial velocity due to the neglect of particle 
inertia. Predictions this close to the injector exit are extremely 
sensitive to initial conditions, however, which may be a factor in 
these results. 

The radial profile of particle mass flux illustrates the 
differences between the three models. In contrast to the nonswirling 
flow, the peak particle mass flux for the swirling flows is not found 
at the centerline of the jet, but is shifted radially outward due to 
centrifugal forces. Calculations for all three models were initiated 
at x/d * 0.5 using the experimentally measured particle mass fluxes 
as initial conditions. As illustrated in figures 4.44 and 4.45, the 
locally homogeneous flow model predicts a shift in the maximum mass 
flux to the center of the jet, because turbulent dispersion of the 
particles is overestimated due to the no-slip assumption. Predictions 
of particle mass flux from both of the separated-f low models are 
similar and show better agreement with measurements. Stochastic 
separated flow model predictions of particle fluctuating velocities 
at x/d = 2 are also illustrated in figures 4.44 and 4.45 and show 
good agreement with measurements. Even though continuous-phase 
fluctuations are assumed to be isotropic, predicted fluctuating 
particle axial velocities are greater than fluctuating radial or 
angular velocities because particles with different axial velocities 
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are transported radially by fluctuating radial and angular 
velocities. Also illustrated in figures 4.44 and 4.45 are particle 
angular velocities. Angular velocities at x/d = 2 appear to be 
better predicted using the locally homogeneous flow model, however, 
this may be the result of inaccurate initial gas-phase angular 
velocities or the assumption of equal initial angular velocities of 
all size groups, as discussed earlier. Measured angular velocities 
also have relatively high uncertainties, as discussed in Appendix A. 

Radial profiles of number-averaged flow properties for the 
swirling, particle-laden jets at x/d = 5 are illustrated in figures 
4.46 and 4.47 for S = 0.16 and 0.3, respectively. As expected, 
increasing the swirl number increases the width of the particle-laden 
jet, similar to the behavior observed for the single-phase jets. The 
gas-phase flow width is slightly overestimated at x/d = 5 for both 
swirling flows. Again, this is probably due to turbulence modulation 
by the particles, as discussed earlier. Also, because the gas-phase 
flow width is overestimated, particle axial velocities are also 
overestimated for the separated-f low models. The stochastic separated 
flow model overestimates particle axial velocities to a greater extent 
than the deterministic separated flow model. The locally homogeneous 
flow model, however, underestimates particle axial velocities, due to 
the neglect of particle inertia. At x/d = 5, particle axial 
velocities, number averaged over all size groups, are nearly equal to 
the continuous phase. A comparison between predictions and 
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measurements of particle mass flux at x/d = 5 shows similar results 
to those observed at x/d = 2. However, differences between the model 
predictions are more pronounced. For both swirling flows, the locally 
homogeneous flow model predicts the maximum particle mass flux at the 
center of the jet, which is clearly not correct. The neglect of 
turbulent dispersion of the particles, illustrated by the 
deterministic separated flow model, causes the particles to be 
confined to a relatively narrow region of the flow. If only mean 
properties of the continuous phase are considered, the particles are 
transported by centrifugal forces to regions where radial velocity is 
small and tend to remain there. Considering turbulent fluctuations of 
the continuous phase gives better predictions of particle mass flux 
at x/d = 5. The predicted maximum particle mass flux is shifted 
radially outward when compared to the measurements. This is again 
probably caused by the overestimation of jet width at this streamwise 
location. Predictions of mean particle angular velocities and 
fluctuating particle velocities are quite good for the S = 0.16 flow, 
see figure 4.46. For the higher swirl number flow, S = 0.3, particle 
fluctuating axial velocities are overestimated, however, predictions 
of the other two fluctuating particle velocities show better agreement 
wi th measurements . 

At the streamwise location of x/d = 10, number-averaged particle 
velocities, averaged over all sizes, are greater than the continuous 
phase for both swirling flows considered during the present study. 
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Radial profiles of flow properties are illustrated in figures 4.48 and 
4.49 for S = 0.16 and 0.3, respectively. The continuous-phase flow 
width is again slightly overestimated, however, predictions show better 
agreement with measurements at x/d = 10 than at locations closer to 
the injector. For the lower swirl number flow, S = 0.16, particle 
axial velocities are well predicted with the stochastic separated flow 
model. Predictions based on the locally homogeneous flow model also 
are quite good, however, this model overestimates the rate of decay of 
axial velocities see figure 4.42; therefore the unnormalized predicted 
velocities are lower than the measurements. At this streamwise 
distance, angular velocities have decayed to negligible values and are 
not reported. As illustrated in figures 4.48 and 4.49, the radial 
distribution of particle mass flux at x/d = 10 has significantly 
shifted from the profiles found at x/d = 5. By x/d = 10, the 
particles have dispersed sufficiently so that the maximum mass flux 
has shifted to the center of the jet for the S = 0.16 flow and is 
clearly headed in that direction for the S = 0.3 flow. Neglecting 
turbulent dispersion of the particles causes the particles to remain 
in a relatively narrow region for both swirling flows. For the 
S = 0.3 flow, no particles whatsoever are predicted for r/x less 
than about 0.1 with the deterministic separated flow model. The 
stochastic separated flow model does not adequately predict this abrupt 
shift in particle mass flux at this streamwise location. Predicted 
values of fluctuating radial and angular number-averaged particle 
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velocities are lower than measurements for both swirl flows. This is 
probably due to the assumption of isotropic velocity fluctuations for 
the eddies in the stochastic model. 

Radial profiles of gas-phase properties at x/d = 20 are 
illustrated in figures 4.50 and 4.51 for both swirl flows. Predictions 
are in better agreement with measurements at this position than closer 
to the injector, although the jet width is still slightly 
overestimated. Predictions at x/d = 20 are not as sensitive to 
initial conditions and the swirl component has almost completely 
decayed. Again, there is little difference between predictions of 
gas-phase properties for the no-slip and separated-f low models. 
Predictions of axial velocities are in good agreement with measurements 
for all three models. However, since the locally homogeneous flow 
model overestimates the rate of decay of axial velocities, predicted 
unnormalized velocities from this model are lower than the 
measurements. Predicted values of (u ' ^) 1 /2 underestimate the 
measurements while (v ,2 ) 1/2 and (w' 2 ) 1/2 are in reasonably good 
agreement with measurements using the stochastic model. Since effects 
of swirl have decayed at this axial location, ignoring the anisotropy 
of the continuous phase is the main reason for this behavior. The 
particle mass flux predictions again highlight the different physical 
assumptions of the three models. Particle mass flux measurements 
indicated that between x/d = 5 and x/d = 10, the maximum mass flux 
shifted to the center of the jet for both swirl cases. Since angular 
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and radial velocities have decayed to relatively small values at this 
distance, and would tend to move particles outward, the only mechanism 
for transport inward is turbulent dispersion. As shown in figures 
4.50 and 4.51, the predicted maximum particle flux from the stochastic 
separated flow model has not completely shifted to the center of the 
jet, however, it is clearly evolving in this direction. In contrast, 
the deterministic separated flow model predicts a very narrow 
distribution with no particles at the center of the jet. Compared to 
the nonswirling case, the locally homogeneous flow model 
underestimates particle dispersion for both swirling flows at x/d = 

20. This behavior is caused by neglecting the angular inertia of the 
particles which tends to transport them radially. 

Measurements and stochastic separated flow model predictions of 
mean axial velocities, for particle diameters of 23, 43, and 63 M m in 
the two swirling flows, are illustrated in figures 4.52 and 4.53. 

Radial profiles of u p are shown at axial locations of x/d = 2, 5, 

10, and 20. At x/d =2, the velocity of the smaller particles is 
larger than the larger particles near the center of the jet. Due to 
their increased inertia, the axial velocities of the larger particles 
decay at a slower rate than smaller particles, and by x/d = 20, 
larger particles are moving at higher velocities than the smaller 
particles at all streamwise locations. Increasing swirl increases the 
variation of velocity with particle size, see figures 4.52 and 4.53. 
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Stochastic separated flow model predictions of mean particle 
axial velocity for each size group show better agreement with 
measurements as distance from the injector exist increases. At 
streamwise distances of x/d = 10 and less, radial profiles of axial 
velocity are overestimated for both swirling flows. As discussed 
earlier, this is probably due to turbulence modulation since the 
continuous-phase jet width was also overestimated for these flows, see 
figures 4.44 to 4.49. At x/d = 5, particle velocities are 
underestimated for both swirling flows. Since predicted centerline 
particle velocities were higher for the deterministic separated flow 
model at x/d = 5 and 10, this appears to be the result of eddy 
specification in the stochastic model. This is especially true for 
the 23 pm particles. Eddies are assumed to travel at the gas-phase 
velocity at their point of origination; therefore, the smaller 
particles, which have shorter relaxation times, tend to remain in a 
particular eddy longer than larger particles. 

Predictions from the stochastic model display the correct trends 
with respect to changes of particle velocities with streamwise 
distance. Predictions at x/d = 20 exhibit good agreement with 
measurements for both swirling flows. 

4.3 Sensi ti vi tv Study 

Based upon the results of the present study, the stochastic 
separated flow model appears to be reasonably successful in treating 
weakly swirling, particle-laden jets. All predictions, however, 


147 


depend on the prescription of initial conditions and the empirical 
correlations for Cp. The sens! ti vi ty of the stochastic separated 
flow model to variations in various parameters is examined in the 
present section. 

Nine key parameters were considered during the sensitivity 
study. The initial turbulence kinetic energy was considered since it 
could not be measured for the particle-laden jets and single-phase jet 
measurements were used instead. Similarly, Cq was computed using 
k Q and a length scale based upon the jet width; therefore, effects 
of uncertainties of e Q were also studied. Since the present flows 
were not monodi sperse , d p was also chosen as a variable in the 
sensitivity study. Standard particle drag empirical relations were 
adopted, thus, Cp was also considered. The eddy length scale, 
which was adopted from Gosman and Ioannides (66), as modified by Shuen 
(41), was also considered. Since initial particle radial and angular 
velocities were number averaged over all size groups, these were 
considered as well. For completeness, initial particle axial 
velocities were also considered. Finally, since gas-phase angular 
velocities could not be measured in the particle-laden jets, and were 
estimated by subtracting particle angular momentum from single-phase 
measurements, initial values of w were also considered during the 
sensitivity study. 

Results of the sensitivity study are summarized in tables 4.1 and 
4.2. The percent change in the computed output variable for a 
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Table 4.1 - Results of the Sensitivity 
Study at x/d = 5 for S = g.3 
Particle-Laden Swirling Jet 


Param- 

eter 

u 

c 

k 

c 

w 

m 

u 

pm 

w 

pm 

b 

G 

m 


u 2 
c 




k 

0 

-0.8 

2.5 

-1.6 

-1.8 

-3.5 

-0 

c 

0 

.3 

-.6 

.6 

.8 

3.2 

-4 

w 

0 

-5.6 

9.5 

19.5 

-10.0 

4.5 

-0 

u 

po 

-0 

-0 

-0.1 

-2.8 

1 .2 

-0 

w 

po 

~0 

.2 

1 . 1 

-9.1 

-5.6 

8 

V 

po 

~0 

-0 

-.2 

-10.4 

-1.3 

-0 

C D 

. 1 

-.1 

-.3 

-.6 

3.6 

-17 

d 

P 

.7 

-2.5 

-1 .9 

46.0 

-6.6 

-8 

i 

L 

e 

-. 1 

.3 

-.6 

3.3 

.3 

8 


d Percent change in output variable for a 
25 percent increase in input variable. 
^Percent change in radial location of G, 








25 percent increase in the input variable are tabulated. The 
sensitivity study was conducted using the stochastic separated flow 
analysis. Baseline predictions were generated using initial conditions 
of the S = 0.3, particle-laden jet, with a single particle size of 39 
mi crons . 

At x/d = 5, see table 4.1, a 25 percent increase in any of the 
input parameters has a relatively small effect on all gas-phase 
predicted properties, except for w . This is not surprising since 
jet structure is very sensitive to swirl number. Increasing “w^ by 
25 percent increases both the turbulence kinetic energy and maximum 
angular velocity. A 25 percent increase in w q causes a 5 percent 
decrease in centerline axial velocity since the width of the jet is 
increased. Gas-phase predicted properties are not greatly influenced 
by initial particle parameters since this flow is relatively dilute. 

Predictions of particle properties at x/d = 5 show a greater 
dependence on initial conditions than the continuous phase. As 
illustrated in table 4.1, increasing the particle diameter by 25 
percent produces the largest effect on predictions of particle 
properties at x/d = 5. The increased particle diameter causes the 
predicted mean particle axial centerline velocity to increase by 
nearly 50 percent due to the increased inertia of the larger 
particles. Predictions of particle properties are less sensitive to 
initial values of particle angular and radial velocities. Increasing 
w" causes the particles to move to larger radial locations where 
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they encounter reduced continuous-phase axial and angular velocities, 

tending to reduce particle velocities. Increasing w also causes 

an increase of the radial location of G... Increasing 7 also 

M po 

decreases particle velocities, however, G., is not affected. 

Variations in the other parameters considered during the sensitivity 
study produced smaller changes in particle properties at x/d = 5, see 
table 4.1 . 

Table 4.2 presents the results of the sensitivity study at 
x/d = 20. For the baseline case considered, all velocities have 
decayed to relatively low values, especially the angular velocities of 
both phases, which are negligible at this streamwise location and are 
not reported. As expected, sensitivity of the predictions to changes 
in initial conditions is reduced at x/d = 20 compared to x/d = 5, 
with the exception of the radial location of G w . At x/d = 20, the 
radial location of G M is shifted toward the center of the jet, so 
that small variations in the predicted location show up as large 
percentage changes. It is evident from table 4.2 that the predictions 
show the largest sensitivity to the particle diameter but the increase 
in particle axial velocity is reduced to about 5 percent at x/d = 20. 
Predictions are less sensitive to the other parameters at x/d = 20. 

For the weakly swirling particle-laden jet considered during the 
sensitivity study, the stochastic separated flow predictions are most 
sensitive to particle size. Predictions are less sensitive to initial 
values of "w, v7 and 7 although these parameters are still 

H r 

important. Predictions appeared relatively insensitive to the other 
parameters considered during the sensitivity study. 



CHAPTER V 


SUMMARY AND CONCLUSIONS 


5. 1 Summary 

The overall objective of the present study was to investigate 
weakly swirling, particle-laden, turbulent jets. Measurements were 
emphasized, however, predictions were used to help interpret the 
measurements and to initiate evaluation of methods to estimate flow 
properties . 

Experiments were initially conducted for three, single-phase jets 
having swirl numbers of 0, 0.19, and 0.33 to provide baseline data for 
the particle-laden jets. Measured flow properties for the 
single-phase jets included: mean and fluctuating axial and angular 

velocities, fluctuating radial velocity, turbulence kinetic energy, and 
Reynolds stress. A two-color LV system was used to measure velocities 
of the single-phase jets. 

Experiments were also conducted for three particle-laden jets 
with swirl numbers of 0, 0.16, and 0.3. The particle size 
distribution had a SMD of 39 pm while a single mass loading ratio of 
0.2 was used. For the continuous phase, mean and fluctuating axial 
velocities were measured with a single-channel phase/Doppler 
anemometer. For the particle phase, mean and fluctuating axial, 
radial, and angular velocities were measured using a two-color LV 
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system. Also for the particle phase, mean and fluctuating axial 
velocities for each particle size were measured using a single-channel 
phase/Doppler anemometer. Isokinetic sampling was used to measure 
particle mass flux distributions. 

For the single-phase flows, governing equations for conservation 
of mass and momentum were solved using the finite-difference code, 
GENMIX (89). The equations were closed using a k-c turbulence model 
which had been calibrated for constant and variable density, 
single-phase, round jets. Two modifications of the k-c model, which 
attempt to account for effects of streamline curvature, were also 
evaluated. 

Three two-phase flow models were evaluated for the particle-laden 
flows. They included: (1) a locally homogeneous flow (LHF) model, 

where interphase transport rates are assumed to be much faster than 
the rate of development of the flow as a whole; (2) a deterministic 
separated flow (DSF) model, which allows for finite interphase 
transport rates (evaluated using the mean properties of the continuous 
phase), but ignores interactions between the particles and turbulent 
fluctuations; and (3) a stochastic separated flow (SSF) model, where 
finite interphase transport rates and interactions of particles with 
turbulent eddies are considered using Monte-Carlo methods. 

The governing equations for the continuous phase were based on 
the Favre-averaged conservation equations written in an Eulerian 
coordinate system. The dispersed phase was treated (for the separated 
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flow models) by solving Lagrangian equat ions-of-mot ion for the 
particles. A modified version of the GENMIX program (89), combined 
with a second-order Runge-Kutta ordinary differential equation solver 
for particle motion, was used to solve the governing equations. 

5.2 Conclusions 

The major observations and conclusions of the present study are 
as follows: 

(1) For the weakly swirling free jets, increasing the swirl 
number increases the rate of mean axial and angular velocity decay 
with streamwise distance. Increasing the swirl number also increases 
the turbulence kinetic energy, the Reynolds stress, and the width of 
the swirl ing jets . 

(2) For the single-phase, weakly swirling jets, a version of the 
k-e turbulence model, which was modified to include effects of 
streamline curvature, showed better agreement with measurements than 
the standard k-e turbulence model. This modification involved 
replacing one of the constants in the dissipation equation with a 
function of the flux Richardson number. 

(3) Predictions using the stochastic separated flow model showed 
reasonable agreement with measurements for the nonswirling, 
particle-laden jets. In general, the locally homogeneous flow model 
overestimated both the rate of particle velocity decay and the rate of 
spread of particles in the nonswirling jet due to the neglect of 
particle inertia. The deterministic separated flow model 
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underestimated particle spreading rates due to the neglect of the 
effect of turbulent fluctuations on particle motion. Only the 
stochastic separated flow model, which accounts for both particle 
inertia and effects of turbulent fluctuations, correctly predicted 
particle spreading rates over the entire flowfield. 

(4) Particle axial fluctuating velocities were generally 
underestimated at far downstream locations. This is probably due to 
the assumption of isotropic velocity fluctuations in the stochastic 
separated flow model, since fluctuating axial velocities are expected 
to be greater than fluctuating radial on angular velocities. 

(5) Near the injector exit, jet widths were overestimated with 
the separated-f low models. This was probably caused by turbulence 
modulation by the particles, which was not considered in the analysis. 

(6) Mean axial velocities for each particle size group were 
reasonably well predicted for the particle-laden jets using the 
stochastic separated flow model. Predictions showed the same trends 
as measurements for particle streamwise axial velocity decay of each 
size group. 
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APPENDIX A 


EXPERIMENTAL UNCERTAINTIES 

A . 1 Bias Erro rs 

A.1.1 Sin gle -Phase Jet Velocit ies 

Bias errors in velocity determinations using LV arise ftom 
several sources, as follows (93): (1) directional ambiguity, due to 

the inability of a stationary or slowly translating fringe pattern to 
provide an indication of the direction in which particles are crossing 
the fringes; (2) directional bias, due to particles crossing the 
measuring volume at a small angle with respect to the plane of the 
fringes so that an insufficient number of fringes are crossed to be 
processed; (3) concentration bias, due to varying particle 
concentrations in the flowfield; (4) velocity bias, due to the fact 
that for a uniformly seeded flowfield, more particles having a higher 
velocity are measured than those having a lower velocity; and (5) 
gradient bias, due to variations in flow velocity in the measuring 
volume. Each of these possible sources of biasing are discussed in 
the following. 

The present measurements in single-phase jets employed frequency 
shifting for both channels to eliminate errors due to directional 
ambiguity. Effective shifting levels were set at each location so 
that maximum negative velocities could be detected. For axial 
velocity measurements, effective shifting levels, of 0.5 and 1 MHz 
were used. For radial and angular velocity measurements, effective 
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frequency shifting levels of 2 and 5 MHz were used. Use ol frequenry 
shifting also minimized directional biasing. 

Concentration bias was not a factor during the present 
experiments since both the jet and ambient surroundings were seeded. 

Velocity bias can be a factor in highly turbulent flows when 
particle averages are used. Various techniques have been suggested 
(93 to 97) to minimize this effect. The approach taken during the 
present measurements was to operate the counters in the total burst 
mode and apply a weighing factor to the measurement based upon the 
measured time of the doppler burst. This approach was suggested by 
Buchave and George (93) and found to give reasonably good results. 

For the present measurements, time averaged values were calculated as 
fol lows : 
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where x . is the measured time of the ith doppler burst. Time 
averaged values of mean and fluctuating, radial and angular velocities 
are calculated in a similar manner. For Reynolds stress, time 
averaged values were calculated as follows: 
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where x . is the smaller of the two measured doppler burst 
uvi 

times from the two counters measuring u and v. 

The final source of biasing considered is gradient biasing. This 
bias can be estimated using the analysis reported by Kreid (98). For 
this analysis, only velocity gradients in the radial direction are 
considered and all seeding particles crossing the measuring volume are 
assumed to be recorded, so that higher-order terms in the Taylor 
series expansion of mean velocity (expanded from the center of the 
measuring volume) can be ignored. The difference between the actual 
velocities, u and u', and the measured values, u^. and u ^ , 
can then be estimated as follows: 
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where 4 is the half-width of the measuring volume in the radial 
direction . 

Maximum gradient biasing occurred at x/d = 0.5 for the flows 
considered here. Using the depth of field limit for measuring volume 
length, the estimated gradient biasing error was less than 0.1 percent 
for mean values and approximately 0.5 percent for fluctuating values. 
If the depth of field limit for probe volume length is not applied, 
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gradient bias errors are larger, however, and at the same location, 
they are estimated to be less than 1 percent for mean values and 
approximately 30 percent for fluctuating values. These are maximum 
values, and are only found in a small radial region at x/d = 0.5: 
the majority of errors are considerably smaller at other locations in 
the flow. 

A.1.2 Particle-Laden Jet Vel oci ti es 

A. 1 ,2 . 1 Pi sp ers ed -p hase velocities - Biasing errors for LV 

measurements of particle velocities are similar to single-phase 
velocity measurements. 

Frequency shifting was used to eliminate errors due to 
directional ambiguity and to minimize errors due to directional bias. 

Velocity bias was not a factor in these measurements because 
number-averaged velocities are desired for comparison with theory. 

For particle velocity measurements, the counters were operated in the 
total burst mode but no residence time weighting factor was included 
so that number averages were obtained. 

Gradient biasing was the greatest potential source of bias error 
for particle velocity measurements. For these measurements, grazing 
encounters with the measuring volume are recorded as valid 
measurements so that the dimensions of the measuring volume were 
increased by the diameter of the largest particles in the flow. 
Equations (A. 4) and (A. 5) were used to estimate gradient bias errors. 
Maximum bias errors are found at x/d = 0.5. For the worst case 
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considered, particle velocities are biased less than 0.1 percent for 
mean values and less than 2 percent for fluctuating values. At other- 
positions in the flow, biasing errors for fluctuating values are 
reduced to less than 1 percent. Therefore, gradient biasing errors 
are small in comparison to other uncertainties of the particle 
velocity measurements. 

A. 1.2. 2 Continuous-phase v elocities - Biasing errors for LV 
measurements of continuous phase velocity are similar to those 
previously discussed. The instrument used to measure continuous phase 
velocities simultaneously measured both particle size and velocity. 
Because frequency shifting was not available for this instrument, only 
mean and fluctuating axial velocities were measured in order to 
minimize directional bias. To further minimize directional bias, 
axial velocity measurements were not attempted in regions of the flow 
near the edge of the jet where negative velocities may be found due to 
the low mean values of axial velocity. 

Velocity bias may influence measurements of continuous-phase 
velocity, especially in highly turbulent regions. The processor 
considered the entire Doppler burst and mean and fluctuating axial 
velocities were calculated from number-averaged velocities for each 
particle size group. As discussed in (93, 94, 96), for flows with 
turbulence intensities above approximately 30 percent, velocity bias 
can be an important source of error in measurements. For the present 
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measurements, regions near the edge of the jet where turbulence 
intensities are quite high were avoided in order to minimize velocity 
bias effects. 

Gradient biasing errors for continuous-phase measurements are 
similar to those previously discussed and were estimated using 
equations (A. 4) and (A. 5). Maximum biasing occurred at x/d = 0.5 
where gradient bias errors were estimated to be on the order of 0.1 
percent for mean values and 1 percent for fluctuating values. It can 
be concluded that gradient bias was negligible for the present 
measurements . 

A. 1 .3. Parti c le Mass Flux Measurement s 

Particle mass flux measurements were obtained by collecting 
particles using an isokinetic sampling probe for a timed interval and 
weighing. In order to insure reasonable sampling times and adequate 
resolution, probes with inside diameters of 2 and 5 mm were used. 
Because of the relatively large diameter of the probe, errors due to 
particle gradients in the radial direction are the major sources of 
bias errors. Similar to the analysis for gradient bias in velocity 
measurements, the gradient bias error in particle mass flux 
measurements can be estimated as follows: 



(A. 6) 
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where d is the inside diameter of the probe. Gradients ate highest 
at x/d = 0.5 where the maximum gradient bias was estimated to be on 
the order of 2 percent. 

A . 2 Uncertainty Estimates 

A.2.1 General Formu lation 

The uncertainty analysis described by Kline and McClintock (99) 
and Moffat (100) was adopted to estimate experimental uncertainties 
for the present measurements. An output variable, R, is consideted, 
which is a function of n measured variables, X., as follows. 

R = R(X 1 ,X 2 ,X 3 ,...X n > <A. 

The X. are subject to uncertainties, Sx., and the resulting 
uncertainty in R, e.g., SR, is to be determined. The SX. and 
6R can be defined as the expected standard deviation of these 
quantities as some percentage of the confidence interval of these 
quantities, e.g., the 95 percent confidence interval is most 
frequently chosen (100) and is used here as well. The relative 
uncertainty can be expressed as follows: 

SR I" 9R ft* A __ 3R_ 

R " 9X ] \ R / + 9X 2 

Equation (A. 8) will be used to estimate relative uncertainties in the 
fol lowing. 

A.2.2 Mean and Turbulent Continuous-Phase Velocities 
A.2.2.1 Single-phase jet measurements - In considering 
uncertainties in continuous-phase veloci ty measurements , positioning 
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errors are ignored since this only influences the position at which 
data points are located and this accuracy has been given previously. 

The mean velocity measurement then depends on the overall calibration 
factor, K, between the averaged electrical signal, E, from the counter 
and velocity. The electrical signal, E, is averaged over all Doppler 
bursts considered. This can be expressed as follows: 

U = KE (A. 9) 

x 

Applying equations (A. 8) to (A. 9) yields: 

3u 
u 

The overall calibration factor, K, was verified using a rotating 
disc. The uncertainty of E was estimated from actual measurements 
at various locations in the flowfield using standard operating 
procedures. The uncertainty is estimated as two times the standard 
deviation found in these measurements <~95 percent confidence 
interval). Substituting the appropriate estimates into equation 
(A. 10) then yields an uncertainty for mean axial velocities of 
5 percent. Uncertainties of mean angular velocity measurements can be 
quite high since values of w decrease to zero at the centerline and 
near the edge of the jet. In addition, peak values of w are smaller 
than TT and w decays rapidly with axial distance. Using equation 
(A. 10) for w at x/d = 10 yields an uncertainty for mean angular 
velocity of 20 percent. Uncertainties in w are smaller closer to 




the injector exit, since w is larger while uncertainties farther 
downstream are larger because w is smaller. 

The uncertainties in velocity fluctuations were also estimated 
for the single-phase flows. For this case e^ represents the 
fluctuating electrical signal from the doppler bursts and fluctuating 


axial velocity can be expressed as follows: 

u' . Ke x (All) 

Applying equation (A. 8) yields 



Equations for 7' and w' are similar. Using representative 
standard deviations in measurements of fluctuating quantities 
throughout the flowfield yields uncertainties for u, V ' , and w 1 of 
less than 5 percent. 

Uncertainties in quadratic quantities, such as k and u'y', are 
larger. For k, we have: 


k - 1 K 2 (e^ , 5 2 ) 

Assuming e’ ~ = ®q» ar, d ^ppl yi ng equation (A. 8) yields. 


(A. 13) 


Sk 

k 



(A. 14) 


Substituting appropriate estimates into equation (A. 14) yields an 


uncertainty estimate of 12 percent for k. 
For Reynolds stress, we have: 
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v 2~ ~ 
u v •- K c e 
x i 


Applying equation (A. 8) yields 


6u ' v 


u v 




1/2 


(A. 15) 


Obviously, similar to w - , in regions where 'u'v’ 1 is small, (near 
the axis and edge of the jet) the uncertainty can be very large. A 
more representative condition is the region near the peak of maximum 
Reynolds stress. Substituting these values into equation (A. 15) 
yields an uncertainty of 14 percent for peak values of 'u'V . 

A.2 .2 .2 Pa r ticle - Laden Jet M easurements - Uncertainties of 

mean and fluctuating axial velocities of the continuous phase for the 
particle-laden flows were calculated from equations (A. 10) and 
(A. 12). Uncertainties are estimated to be less than 4 percent for XL 
and less than 5 percent for tf ' . These are identical to the estimated 
uncertainties for the single-phase flows. 

A.2.3. Oi spersed-Phase Velocities 

Uncertainties of measured mean and fluctuating particle 
velocities were calculated from equations (A. 10) and (A. 12). The 
estimated uncertainty of u\ was less than 5 percent. The 

estimated uncertainty in V was 10 percent at x/d = 0.5, which is 

r 

larqer than XT because 7 is smaller. The estimated 

uncertainty in w at x/d = 0.5 is 6 percent near the maximum 

value of w , however, for the same reasons as previously discussed 
P 
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for other measurements of angular velocity, uncertainties in w p 
can be nearly 20 percent since w p varies considerably throughout 
the flowfield. For fluctuating particle velocities, uncertainties are 
estimated to be less than 5 percent. 

A.2.4. Particle Mass Flux Measurement s 

Particle mass flux was determined from the following equation: 


G = 


4 ^ 

* d 2 


(A. 16) 


Applying equation (A. 8) to equation (A. 16) yields: 



(A. 17) 


Substituting estimated uncertainties of probe diameter and mass 
flowrate into equation (A. 17) yields an uncertainty estimate of less 
than 8 percent for particle mass-flux measurements. 
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APPENDIX B 


DESCRIPTION OF THE GENMIX ALGORITHM 
A general purpose computer code, GENMIX, described by Spalding 
(89) for boundary-layer flows, was used to solve the governing 
equations . This code is very convenient for the calculation of free 
jets since it utilizes a dimensionless stream function formulation 
that automatically expands the cross-stream grid width as the 
calculation proceeds downstream. The dimensionless stream condition 
i s defined as fol lows : 

(♦ - t i) 


(*• • *i) 


(B. 1 ) 


where tj and ^ are the values of the stream function at the inner 

and external boundaries of the flow. 

The stream function is defind as: 


and 


8x 


8«I> 


rpv 


( B . 2 ) 


8r 


= + rpu 


( B . 3) 


The use of the stream function automatically satisfies the 
conservation of mass equation. At the inner and external boundaries 
the stream function can be expressed as: 

8HG 


I 


8x 


- rjBij 


( B . 4) 
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3x 


r E m E 


( B . 5 ) 


where irij" and m E " are the mass transfer rates across the inner 
and external boundaries. 

Using these transformations, equation (3.2) can be transformed to 
the following general form: 

„ / \ 

( B . 6) 


3<£ 

3x 


ll * A 4 


3< P _ 3 - (a - 3 - 

'1 T n 2“" 3u = 3w \ 3 3, 

where <p represents u, 7, k and c. The other terms in equation 

(B. 6) are defined in table B.l. For <p = rw, equation (B.l) can be 

rearranged to eliminate the source term and can be written as follows: 


(B . 7 ) 


rgW 

3(rw) /A . \ 3( rw) 3 . 2( _r_ 

IT + (A 1 + A 2 u) = 3u A 3 r \3w 

where A,, A 9 and A^ are also defined in table B.l. This set 

of equations was integrated from an upstream initial condition to a 

value of x equal to 30 diameters, for values of w ranging from 

zero to one . 

The present calculations were performed using 33 cross-stream 
nodes. The cross stream grid spacing included the entire flow width. 
The forward marching step was limited so that the quantity of fluid 
entrainment during the mixing step is a certain fraction of the total 
fluid in the flow to that position in the flow. 

(*E ‘ t l) 


(B. 8) 


r i m i - r E m E 
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Tdble B.l. - Definition of Terms in the 
Generalized GENMIX Program 


A zrm / ( ^ - 40 

1 LIE I 


A ={rm"-rm n )/{4' - 40 

2 EE LIE I 


2 — 2 
A = r puu / (t r - 40 o 

j t t 1 4> 


4 > 


8 x 

pu 


f 


k 
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In the present calculations, the forward step size was limited so that 
the ratio of entrainment to the total fluid in the Mow was 0.05 tor 
the nonswirling flows and 0.02 for the swirling flows. The ratio was 
set as low as 0.005 without significant difference in the results. 
Forward step size was also constrained to be less than 5 percent of 
the current flow width for the nonswirling flows and 2 percent for the 
flows wi th swirl. 
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APPENDIX C 
SUMMARY OF DATA 


C. 1 Single-Phase Jets 

C.1.1 Single-Phase Jet, S = 0 

Table C.l - Single-Phasg Jet at 
x/d = 0.5, S = 0 


r 

X 

_y 

u 

c 




u 1 v 1 2 

xio 

u 

c 

0 

1.0 

mmm 



0.008 

. 121 

.997 


.039 

.026 

.022 

.279 

.972 

.056 

.049 


.039 

.321 

.962 

.058 


.033 

.084 

.479 

.917 

.069 

.043 

.042 

.097 

.521 

.905 

.074 

.045 

.048 


.679 

.848 

.077 


.052 


.721 

.827 

.078 

.052 

.058 


.879 

.721 

! . 1 16 

.067 

. 1 12 

.337 

.923 

.663 

. 127 


. 144 

.480 

1 .017 

.464 

.159 


.227 

.874 

1 .079 

.256 

. 141 

. 118 

.239 

1.027 

1 . 120 

.193 

.115 

.113 

. 193 


1.175 

. 106 

.086 


. 144 

.468 

1 . 180 

.099 

.085 

.101 

.137 

.478 


a u c = 14.86 m/s. 


Table C.2 - Single-Phase Jet at x/d = 15, S = 0 a 


r 

X 

- 

- 2 1/2 
( u ' ) 

- 2 1/2 
(v 1 ) 

(w’ 2 ) 1/2 

k _ 1 

u 1 V 1 2 

=~y~ xio 
u c 
c 

u 

c 

u 

c 

u 

c 

u 

c 

- 7 x 1 0 
U c 

c 

0 

1.0 

0.209 

0. 165 

0.165 

0.489 

0.021 

.02 

.964 

.212 

. 166 

.169 

.504 

.523 

.04 

.843 

.214 

. 165 

. 177 

.521 

1 .068 

.06 

.684 

.212 

. 162 

.177 

.513 

1.283 

.08 

.529 

.197 

. 169 

. 182 

.501 

1 .310 

. 10 

.374 

. 164 

. 148 

. 159 

.370 

.971 

.12 

.257 

. 142 

. 130 

.134 

.276 

.676 

. 14 

.172 

.112 

. 120 

.137 

.229 

.515 

.16 

. 106 

.086 i 

. 108 

.093 

.138 

.338 

. 18 

.061 

.059 ! 


.067 

.083 

.238 


a u c = 6.77 m/s. 
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Table C. 3 


Single-Phase Jet at x/d = 30, S - 0 a 


r 

X 

- 

- 2 1/2 
L u.LJ_ 

- 2 1/2 

— 

U 

c 

2 1/2 
(w 1 ) 

u 

c 

-- ? k\0 ] 
u 

c 

' in 2 

- i x 1 U 
u c 
c 

u 

c 

u 

c 

0 

.02 
.04 
.06 
.08 
. 10 
.12 
. 14 
.16 

1 .0 
.964 
.863 
.741 
.582 
.439 
.314 
.222 
.120 

0.244 
.254 
.247 
.241 
.231 
.204 
. 179 
. 141 
.107 

0.210 

.232 

.217 

.228 

.223 

.206 

.188 

.168 

.146 

0.210 
.224 
.224 
.221 
.243 
.206 
. 168 
. 147 
.122 

0.737 
.841 
. 790 
.795 
.813 
.633 
.478 
.348 
.239 

0.012 
.937 
1.004 
1 . 380 
1 . 695 
1.520 
1.340 
.941 
.448 


a u c = 3.03 m/s. 


C.1.2 Single-Phase Jet, S = 0.19 


Table C.4 - Single-Phase Jet at x/d = 0.5, S - 0.19 a 


r 

X 

u 




w 

w 

m 

-i- xl0 ' 


- 2 1/2 
(u< 1 

- 2 1/2 
(V 1 1 

u 

c 

2 1/2 
(W ) 

U 'V.'. ,n 2 

- n X 10 
u 

c 

u 

c 

u 

c 

u 

c 

U 

c 

o 

1 .0 

0.082 

0.149 

0.149 

0.046 

0.257 

0.018 

) 

991 

.088 

. 146 

.143 

.237 

.247 

.029 

2 

.983 

.087 

.153 

.124 

.500 

.232 

.056 

3 

.976 

.087 

. 148 

.109 

.632 

.208 

.087 

.4 

.965 

,087 

.139 

.092 

.742 

.176 

.121 

. 5 

.951 

.087 

.137 

.088 

.814 

.171 

. 139 

.6 

.937 ! 

.088 

. 114 

.091 

.856 

. 145 

. 177 

. 7 

.91 1 

.096 

. 109 

.097 

.907 

.152 

. 335 

.8 

.878 

. 119 

.107 

. 100 

.948 

. 178 

.500 

9 

.81 1 

.158 

.117 

. 108 

1.0 

.253 

. 875 

1 . 0 

.672 

.207 

.139 

.135 

.876 

.401 

1 . 430 

1 . 1 

.476 

.221 

.171 

. 148 

.638 

.499 

1 .850 

1 .2 

.234 

. 167 

. 154 

.136 

.393 

.351 

1 .230 

1 .3 

.109 

. 102 

.125 

.092 

.136 

. 172 

. 500 


a u c = 12.9a m/s, W m = 3.158 m/s. 
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Table C . 7 - Single-Phase Jet at x/d = 10, S = 0.19 a 




- 2 1/2 

- 2 1/2 

2 1/2 


4 2> „' 

— 


r 

X 





Lx-'—i 

Lwl 1 

w 

U 1 Y 

— ■ i 

, . c. 

xlO" 


u 

u 

U 

U 

w 

U 

u 



c 

c 

c 

c 

m 

<: 

c 


0 

1.0 

0 . 265 

0.240 

0.240 

0.090 

0.921 

0 

K>0 

.027 

.982 

.276 

.249 

.242 

.092 

.984 


.81 

.053 

.926 

.284 

.243 

.241 

.249 

. <190 

1 

12 

.080 

.818 

.274 

.248 

.257 

.552 

1 .01 1 

1 

.48 

.107 

.716 

.268 

.237 

.247 

.514 

.947 

1 

. 73 

. 134 

.609 

.261 

.221 

.221 

.730 

.831 

2 

.21 

. 160 

.489 

.240 

.212 

.218 

.735 

.752 

1 

.94 

. 187 

.382 ! 

.211 

. 199 

. 191 

1 .00 

.604 

1 

.50 

.214 

.289 

. 183 

. 184 

.170 

.849 

.482 

1 

32 

.241 

.210 

.161 

.155 

.160 

.415 

.377 

1 

.03 

.267 

. 155 

.138 

. 149 

. 120 

. 121 

.279 

, 

.07 


a u c = 4.751 m/s, w m = 0.422 m/s. 


Table C. 8 - Single-Phase Jet at x/d = 20, S = 0.19 a 










- 

- 2 1/2 
(u‘ 1 

- 2 1/2 
(v* ) 

, 2 1/2 
(w‘ ) 

k in 1 

u 1 v 1 2 

xlO 
u 

c 


u 

c 

u 

c 

u 

c 

u 

c 

- 7 x 1 0 
u 

c 

0 

1 .0 

0.274 

0.219 

0.219 

0.866 

0.52 

.02 

.969 

.279 

.226 

.222 

.890 

1 .00 

.04 

.926 

.277 

.220 

.227 

.883 

1.51 

.06 

.841 

.270 

.213 

.232 

.860 

2.24 

.08 

.761 

.269 

.229 

.229 

.887 

2.05 

.10 

.664 

.278 

.218 

.215 

.853 

2. 14 

.12 

.540 

.244 

.199 

. 193 

.681 

1 . 74 

.14 

.449 ! 

.212 

.21 1 

.213 

.672 

1.34 

.16 

.435 

.200 

.175 

. 182 

.518 

1.62 

.18 

.31 1 

. 191 

.178 

.155 

.459 

1.25 

.20 

. 197 

.158 

.155 


— 

.99 


a u c = 2.72 m/s. 
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19 - Phase/Doppler 
se Veloci ties at 


= 30, 

s = o a 


,ll/2 

u_ 

1 

u 

U 

c 

C 

1 .0 

0.197 

1 .0 

.201 

.973 

.218 

.892 

.215 

.824 

.214 

.770 

.207 

.608 

. 181 

.577 

. 185 

.541 

. 184 

.460 

.153 

.378 

. 142 

.338 

. 134 

.324 

. 1 19 






C.3.2 Particle-Laden Jet, Particle Velocities, S 
Table C.30 - Phase/Doppler Particle 
Velocities at x/d = 0.5, 

S = 0 



r*T#i 








— 

7 1/2 


K 

Average particle size 

u 

p 

<u' 

P 

> 



jjim 

m/s 

m/s 

0 


52.2 

12. 

249 

■ 

1 


201 


12. 

259 




.401 


12. 

136 

; K 

mufm 


.602 


11. 

741 

1 

.044 


.802 


11 

064 

1 

.117 


.936 


10. 

359 

1 

.230 

1 

.069 


8. 

419 

1 

.434 

1 

.136 


4. 

521 

1 

.897 

0 


61.7 

12. 

394 

0 

.989 


.201 


12. 

330 

1 

.050 


.401 


12. 

221 


.998 


.602 


11 

950 


.942 


.802 


11 

497 

1 

.112 


.936 


10 

533 

1 

.374 

1 

.069 


9 

122 

1 

.290 

1 

.136 


4, 

.450 

1 

.721 

0 


71.3 

12, 

.040 

m 

.160 


.201 


12. 

.092 

■ 

.083 


.401 


12. 

.081 

■ 

.976 


.602 


11. 

.821 


.020 


.802 


11. 

.281 

l 

. 142 


.936 


10, 

.721 

i 

.260 

1 

.069 


9, 

.526 

l 

.210 

1 

.136 


1 , 

.663 


.076 

















Table C.31 - Phase/Doppler Par 
Velocities at x/d = 5, 









Table C.31 - C 


r 

X 

Average particle size 

li/n 

u 

P 

m/s 

1/2 

(u’ 2 ) 

P 

m/s 

0 

52.2 

12.687 

0.727 

.027 


12.470 

. 744 

.053 


11.660 

.957 

.080 


10.393 

1 .206 

. 107 


8.661 

1.339 

.134 


7.046 

1 .344 

■ 160 


5.370 

1.267 

. 187 


3.728 

1.108 

.214 




2.218 

.819 

1 

61.7 

12.660 

0.743 

m 


12.460 




11.849 

.881 

HREfSl 


10.608 

1 .392 

.107 


9.061 

1 .246 

.134 


7.696 

1 .216 

. 160 


6.099 

1.199 

.187 


4.286 

1.145 

.214 


2.613 

.709 

0 

71.3 

12.613 

0.713 

.027 


12.419 

.755 

.053 


12.016 

.777 

.080 


10.855 

1.192 

.107 


9.296 

■ 

.134 

i 

8.112 


.160 


6.738 

1 .088 

. 187 


4.904 

1 .212 

.214 


2.637 

.749 




























































Table C.35 - Phase/Doppler 
Gas-Phase Velocities at 
x/d = 2, S = 0. 1 6 a 



Table C. 36 - Phase/Doppler 
Gas-Phase Velocities at 
x/d = 5, S = 0. 1 6 a 




2 1/2 

r 

u_ 

X 

— 

— 


u 

u 


c 

c 

0 

1.0 

0.154 

.054 

.877 

. 157 

.107 

.605 

. 185 

.134 

.451 

. 150 

.160 

.328 

.125 

.187 

.226 

.098 

.214 

. 123 

.065 

.241 

.061 

.032 


a u c = 9.75 m/s. 
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C.3.4 Par ticle-Laden Jet, Particle Velocities, S 
Table C.39 - Phase/Doppler Particle 
Veloci ties at x/d =0.5, 






































Table C.41 - Phase/Doppler Particle 
Velocities at x/d = 5, 

S = 0.16 


r 

X 

Average particle size 
nm 

u 

P 

m/s 

1/2 

2 

(u* ) 

P 

m/s 

0 

14.3 

9.936 

1.164 

.053 


8.837 

1.546 

. 107 


6.788 

1.639 

.134 


5.252 

1 .491 

. 160 


3.981 

1.471 

. 187 


2.582 

1 .076 

.214 


1 .843 

.882 

.241 


1.070 

.677 

.267 


.813 

.544 

0 

23.7 

9.700 

1 . 143 

.053 


8.682 

1 .423 

.107 


6.785 

1.580 

.134 


5.280 

1 .493 

.160 


4.070 

1 .437 

.187 


2.673 

1 .094 

.214 


2.002 

1 . 120 

.241 


1 .084 

.644 

.267 


.790 

.510 

0 

33.3 

9.531 

1 . 105 

.053 


8.632 

1.368 

. 107 


6.904 

1.566 

.134 


5.510 

1.478 

. 160 


4.377 

1 .473 

. 187 


3.038 

1 . 141 

.214 


2.308 

1.071 

.241 


1 . 798 

.728 

.267 


.999 

.592 

0 

42.7 

9.592 

1 .088 

.053 


8.693 

1 . 339 

. 107 


7.178 

1.539 

.134 


5.818 

1.500 

.160 


4.692 

1.474 

. 187 


3.330 

1.171 

.214 


2.585 

1.104 

.241 


1.573 

.923 

.267 


1 . 159 

.657 














































Table C.43 - Cont. 


r 

X 

Average particle size 
*xm 

u 

P 

m/s 

1/2 

2 

(u' ) 

P 

m/s 

0 

52.2 * 

3.931 

0.829 

.02 


3.688 

.845 

.04 


3.448 

.819 

.06 


3.026 

.840 

.08 


2.678 

.816 

. 10 


2.150 

.663 

. 12 


1.781 

.679 

. 14 


1 .477 

.645 

. 16 


1.108 

.466 

.18 


.829 

.359 

0 

61.7 

3.990 

0.826 

.02 


3.724 

.821 

.04 


3.473 

.825 

.06 


3. 1 10 

.849 

.08 


2.703 

.814 

.10 


2.181 

.666 

. 12 


1 .804 

.665 

.14 


1.535 

.640 

.16 


1.161 

.460 

.18 


.854 

.348 

0 

71.3 

3,980 

0.847 

.02 


3.751 

,841 

.04 


3.519 

.820 

.06 


3.116 

.829 

.08 


2.801 

.827 

.10 


2.214 

.650 

. 12 


1.906 

.657 

. 14 


1.618 

.630 

.16 


1.242 

.451 

.18 


.924 

.348 
















C.3.5 Particle-Laden Jet, Gas-Phase Velocities, $ = 0.3 
Table C.44 - Phase/Doppler 
Gas-Phase Velocities at 
x/d = 0.5, S = 0.3 d 


r 

X 

u_ 

u 

c 

, .2.1/2 
(ul 1 

U 

c 

0 

1 .0 

0.074 

.201 

.996 

.081 

.401 

1.017 

.075 

.602 

.992 

.084 

.802 

.975 

.097 

.936 

.778 

.168 

1 .00 

.611 

.205 

1.069 

.310 

.111 

1.136 

.209 

.093 

1.20 

. 142 

.077 


a u c = 1 1 .95 m/s. 


Table C.45 - Phase/Doppler 
Gas-Phase Velocities at 
x/d = 2, S = 0.3 d 






_ 

, ,21/2 

r 


(u 1 ) 

X 

— 

— 


u 

u 


c 

c 

0 

1 .0 

0.096 

.067 

1 .04 

.079 

. 134 

1 .00 

.119 

.201 

.888 

.181 

.267 ! 

.638 

.177 

.334 

.326 

.131 

.368 

.170 

.075 

.401 

. 112 

.062 


a u c = 11.20 m/s . 
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Table C.46 - Phase/Doppler 
Gas-Phase Velocities at 
x/d = 5, S = 0 . 3 a 


r 

X 

u 

C 

•2 1/2 
(y_J 

U 

C 

0 

1 .0 

0.168 

.054 

.909 

.209 

.107 

.733 

.196 

.160 

.476 

. 162 

.187 

.230 

.125 

.214 

.214 

.097 

.241 

.128 

.083 

.267 

.080 

.045 


a u c = 9.35 m/s. 


Table C.47 - Phase/Doppler 


Gas-Phase Velocities at 



A u r =5.75 m/s . 






















r 

X 

Average particle size 

jam 

u 

P 

m/s 

— 1/2 
2 

( u 1 ) 

P 

m/ s 

0 

42.7 

8.825 

1.421 

.201 


8.930 

1.338 

.401 


8.849 

1.175 

.602 


8.512 

1 .071 

.802 


8.292 

1 .064 

.936 


7.868 

1.258 

1.0 


7.390 

1 .435 

1.069 


6.473 

1.544 

1.136 


6.173 


1 .203 


5.721 

HBifl 

1.270 


5.150 

Mmm 

1.337 


4.081 

■EH 

0 

52.2 

8.265 

1.406 

.201 


8.406 

1.309 

.401 


8.394 

1 .086 

.602 


8.221 

.993 

.802 


8. 149 

.984 

.936 


7.854 

1.075 

1.00 


7.462 

1.207 

1.069 


6.917 

1.318 

1.136 


6.737 

1.303 

1.203 


6.436 

1 .233 

1.270 


6.035 

1 .225 

1.337 


5.286 

1.394 

0 

61.7 

7.703 

1.306 

.201 


7.802 

1.238 

.401 


7.783 

1.096 

.602 


7.877 

.948 

.802 


7.968 

.942 

.936 


7.776 

.947 

1.00 


7.507 

1.061 

1.069 


7.057 

1.177 

1.136 


6.968 

1.134 

1.203 


6.724 

1.370 

1.270 


6.521 

1.035 

1.337 


5.942 

1.137 

0 

71 .3 

7.648 

1.329 

.201 


7.714 

1.253 

.401 


7.784 

1.043 

.602 


7.692 

.880 

.802 


7.769 

.900 

.936 


7.481 

.898 

1.00 


7.181 

1.028 

1.069 


6.900 

1.119 

1.136 


6.892 

1 . 140 

1.203 


6.676 

1.117 

1.270 


6.550 

.995 


1.337 


6.057 


1.135 
































r 

X 

Average particle size 

pm 

u 

P 

m/s 

0 

52.2 

9. 147 

.067 


9.396 

.134 


9.463 

.201 


9.211 

.267 


8.454 

.334 


6.947 

.368 


5.358 

.401 


4.300 

0 

61.7 

8.947 

.067 


9.163 

.134 


9.091 

.201 


8.824 

.267 


8.380 

.334 


7.434 

.368 


5.974 

.401 


5.062 

0 

71.3 

8.743 

.067 


8.969 

.134 


8.984 

.201 


8.801 

.267 


8.362 

.334 


7.424 

.368 


6. 149 

.401 


5.325 
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— 1/2 

r 


- 

2 

X 

Average particle size 

u 

(u* ) 



P 

P 


tim 

m/s 

m/s 

0 

52.2 

3.573 

0.773 

.020 


3.496 

.760 

.040 


3.240 

.773 

.060 


2.935 

.776 

.080 


2.502 

.801 

.100 


2.170 

.713 

. 120 


1.771 

.662 

.140 


1 .475 

.628 

.160 


1 .061 

.478 

. 180 


.800 

.391 

0 

61.7 

3.526 

0.796 

.020 


3.523 

.771 

.040 


3.231 

.775 

.060 


2.955 

.788 

.080 


2.560 

.797 

.100 


2.234 

.700 

. 120 


1.808 

.665 

. 140 


1.526 

.631 

.160 


1.090 

.481 

. 180 


.812 

.388 

0 

71.3 

3.629 

.747 

.020 


3.599 

.732 

.040 


3.275 

.776 

.060 


2.971 

.765 

.080 


2.573 

.780 

.100 


2.285 

.699 

.120 


1.911 

.647 

. 140 


1.592 

.623 

.160 


1.192 

.469 
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.382 
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